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ABSTRACT 

One  problem  encountered  in  the  use  of  a  digital 
computer  for  the  purpose  of  simulation  is  the  ordering 
of  simultaneous  events.  This  ordering  is  relative  to 
the  accuracy  of  the  simulation  if  interaction  between 
simultaneous  events  occurs.  In  this  thesis  the  problem 
of  constructing  an  ordering  of  behavior  algorithms, 
called  a  simulation  sequence,  is  discussed.  This  problem 
is  modeled  by  directed  graphs.  An  output  dependency 
relation  and  an  algorithm  dependency  relation  are  defined. 
The  maximal  connected  subset  partition  of  the  algorithm 
relation  graph  is  shown  to  partition  the  simulation  se¬ 
quence  construction  problem.  An  algorithm  splitting 
technique  and  a  selective  search  routine  are  presented. 
These  techniques  are  combined  in  a  general  method  for 
constructing  simulation  sequences. 


iii 


ACKNOWLEDGMENT 


1 

i 


The  author  wishes  to  express  his  debt  of 
gratitude  to  Professor  Herbert  Freeman  for  his  guidance 
and  encouragement  in  the  research  presented,  and  to  Dr. 
Alvy  Ray  Smith  III  for  his  corrections  of  an  initial 
draft  of  this  thesis. 

The  research  reported  here  was  supported  in  part 
by  the  Air  Force  Office  of  Scientific  Research,  Air 
Force  Systems  Command,  USAF,  under  grant  ^p-AF0SR-70-l854 
and  by  the  Bell  Telephone  Laboratories,  Inc.  through  its 
doctoral  support  program. 


iv 


TABLE  OF  CONTENTS 

Page 

List  of  Simulation  Sequence  Properties  vi 

List  of  Illustrations  vii 

List  of  Symbols  x 

I .  INTRODUCTION  2 

1.1  Statement  of  the  Problem  1 

1 . 2  Literature  Survey  2 

1.3  Approach  to  Solution  of  Problem  8 

II.  SYSTEM  MODEL  n 

2.1  Process  Definition  21 

2.2  Network  Definition  24 

2.3  Example  Problem  27 

III.  SIMULATION  SEQUENCE  PROPERTIES  35 

3.1  Basic  Properties  35 

3.2  Structural  Properties  38 

IV.  GRAPH-LEVELING  lj6 

4.1  Level  Sets  llg 

4.2  Application  of  Level  Sets  5I1 

V.  ALGORITHM  ORDERING  61 

5.1  Algorithm  Graph  g2 

5.2  Algorithm  Graph  Partitioning  67 

VI.  ALGORITHM  SPLITTING  74 

6.1  Output  Partitions  75 

6.2  Compatible  Output  Sets  77 

6.3  Maximal  Algorithm  Subset  80 

6.4  Refinement  of  Maximal  Algorithm  Subsets  85 


'  ••  • ‘rTT-Tr*^  sw- 


V 


Page 

VII. 

FORMATION  OF  SIMULATION  SEQUENCES 

92 

7.1 

Modified  Eligible  Set  Method 

92 

7.2 

General  Algorithm 

103 

7.3 

Illustrative  Example 

104 

VIII. 

EXTENTIONS 

119 

8.1 

Single  Algorithm  Characterization 

119 

8.2 

Weighted  Algorithms 

121 

8.3 

Random  Variables 

123 

8.4 

Time  Varying  Network  Structure 

124 

8.5 

Time  Step  Minimization 

128 

IX., 

SUMMARY  AND  CONCLUSIONS 

130 

9.1 

General 

130 

9.2 

Contributions  of  this  Research 

131 

9.3 

Future  Extensions 

132 

X. 

REFERENCES 

133 

APPENDIX 

135 

VITA 

137 

iB"l'.!  '.ilUJSV1,1  ^ 


:^*g?gTgfe»^^^tcyQi?8yBywg^^^*£^igga^fgaE«»<»  «•• 


vi 


LIST  OP  SIMULATION  SEQUENCE  PROPERTIES 


Number 

Property 

Page 

1. 

Covering 

35 

2. 

Minimum  Length 

35 

3. 

Normal  Form 

36 

4. 

”1"  Simulation  State 

37 

5. 

Path  Ordering 

43 

6. 

Existence 

44 

7. 

Minimum  Length 

50 

8: 

Substitution 

55 

9. 

Substitution 

55 

No. 
1 
2 
-  3 

4 

5 

6a 

6b 

7 

8 

9 

10 

11 


vii 

LIST  OP  ILLUSTRATIONS 


System  Model 
Process  Model 

Process  for  which  no  Simulation 
Sequence  Exists 
Illustrative  Example  Network 
Network  Output  Graph  for  Network 
Displayed  in  Figure  A 
Algorithm  Graph  for  Network 
Displayed  in  Figure  4 
Network  Output  Graph  Redrawn  to 
Illustrate  Condensation 
Network  Output  Graph  with  Acyclic 
Algorithm  Dependency  Relation 
Network  Output  Graph  with  Cyclic 
Algorithm  Dependency  Relation 
Maximal  Connected  Subset  Relation 
for  Network  of  Figure  8 
Network  Output  Subgraphs 
Output  Algorithm  Graph  for 
Algorithm  of  Figure  6 


Page 

12 

17 

22 

28 

42 

64 

64 

66 

68 

71 

73 


79 


*3~  -»*• 


viii 


No .  Pag 

12  Output  Algorithm  Graph  82 

13  Network  Output  Graph  87 

14  Network  Output  Graph  with  Maximal 

Subsets  Substituted  for  F^  90 

-15  Flow  Chart  for  the  Selective 

Search  Method  95 


16  Network  Output  Graph  of  Figure  6b 

with  Maximal  Subsets  Substituted 


for  F2  98 

17  Network  Output  Graph  105 

18  Algorithm  Definition  106 

19  Output  Dependency  Relation  Matrix  107 

20  Total  Output  Dependency  Relation 

Matrix  109 

21  Maximum  Algorithm  Definition  111 

22  Network  Output  Graph  with  Maximal 

Algorithm  Subsets  112 

23  Maximal  Subset  Algorithm  Graph 

Matrix  113 

24  Maximal  Connected  Algorithm  Subsets  114 

25  Maximal  Connected  Algorithm  Subset 

Dependency  Relation  Matrix  115 


a-i?'  — ■*  iirf 


>«***SY 


ix 


No. 

Page 

26 

Converse  Level  Partition  of 

the  Maximal  Connected  Algorithm 

Subset  Graph 

116 

27 

Modified  Network  Output  Graph 

122 

28 

Network  Output  Graphs  for  Three 

Network  Configurations 

126 

29 

General  Network  Output  Graph 

127 

=S5JP 


TypCOpW^JMSWES!??^^  '■"~r-  ' 


X 


LIST  OF  SYMBOLS 


Symbol 


a 

a. 


.  a 


M 


a 

a 


R 

N 


$ 

A 

A. 

B 


%  - 

®y  - 


Reference  Page 


algorithm  sequence  16 

output  algorithm  sequence  36 

state  variable  algorithm  sequence  36 

algorithm  sequence  of  length  zero  16 

set  of  input  variables  11 

set  of  network  input  variables  24 

set  of  input  variables  for  process  p  lil 

union  of  all  process  -  input  variable  sets  25 

Input  variable  14 

process  p  input  variable  14 

network  input  variable  24 

empty  set  12 

algorithm  partition  12 

output  algorithm  partition  62 

family  of  maximal  -  connected  algorithm 
subsets  69 

maximal  algorithm  subset  84 

.done  algorithm  set  94 

done  output  set  94 

eligible  algorithm  set  93 

maximal  algorithm  set  94 

eligible  output  set  93 


E 


xi 

Symbol 

Reference  Page 

F 

-  output  algorithm 

15 

f 

-  output  variable  function 

15 

G 

-  state-variable  algorithm 

15 

g 

-  state-variable  function 

15 

H 

-  algorithm 

11 

k 

-  memory-free  function 

14 

K 

-  connected  set  of  algorithms 

67 

L 

-  forward  level  partition 

46 

,.  c 

JLi 

-  backward  level  partition 

52 

£ 

-  number  of  level  sets 

47 

Ky) 

-  level  of  node  y 

48 

M 

-  set  of  state  variables 

11 

m 

-  state  variable 

14 

N 

-  network 

24 

P 

-  set  of  processes 

24 

P 

-  process 

12 

p  ‘vV 

-  path  from  y  to  y 

41 

P(y) 

-  predecessor  set  for  output  y 

86 

R 

-  process  relation 

11 

ra 

-  process  state-variable  dependency 

relation 

15 

rb 

-  process  input  state-variable 

dependency  relation 

15 

R, 


process  state-variable  output 
dependency  relation 


15 


f J*  Vm 


wff’W’gi mmrrr*!***  p-'~ 


Symbol 


Reference  Page 


“  process  Input-output  dependency  relation 
RN  -  network  relation 
Ry  -  output  dependency  relation 

Ry  —  transitive  closure  of  Ry 
RC 

«Y  “  converse  of  relation  Ry 

Rp  ~  algorithm  dependency  relation 

Rk  -  maximal-connected  subset  relation 

S  -  simulation  state  vector 

S(y)  -  successor  set  of  output  y 

si  -  simulation  state  of  variable  i 

V  -  value  vector 

v1  -  value  of  variable  i 

Y  -  set  of  output  variables 

Yp  -  union  of  all  process  -  output  variable  sets 
Yn  ~  set  of  network  output  variables 
y  -  output  variable 


15 

24 

39 

39 

51 

61 

69 

20 

86 

20 

14 

14 

11 


25 


24 

14 


t 


I .  INTRODUCTION 

In  system  design  and  analysis,  simulation  is  fre¬ 
quently  employed  to  s'.udy  systems  that  do  not  lend  them¬ 
selves  to  simple  analytical  methods.  In  circumstances  where 
experimentation  with  the  actual  system  is  impractical,  a 
digital  computer  can  be  used  to  simulate  the  operation  of 
the  system.  Due  to  the  nature  of  the  digital  computer 
several  problems  arise.  Since  a  digital  computer  operates 
in  discrete  time  steps,  the  continuous  time  partitions  of  a 
system  must  be  approximated  by  a  discrete  time  system. 
Furthermore,  since  a  digital  computer  performs  a  single 
operation  at  a  time,  concurrent  events  in  a  system  cannot  be 
simulated  in  parallel.  The  ordering  of  simultaneous  events 
is  the  subject  of  this  dissertation. 

1.1  Statement  of  the  Problem 


The  system  to  be  simulated  is  assumed  to  be  charac¬ 
terized  by  a  collection  (network)  of  subsystems  (processes). 
Associated  with  each  process  is  a  set  of  algorithms  specifying 
the  behavior  of  the  process.  Since  a  Mealy  type  model  (out¬ 
put  is  a  function  of  input  and  state  variables)  is  employed 
for  each  process,  any  change  in  process  status  may  propagate 
instantaneously  throughout  the  network.  Since  the  processes 
interact  the  order  of  execution  of  process  algorithms  in 
updating  the  simulation  time  becomes  relevant  to  the  accurate 
simulation  of  the  system. 


A  simulation  sequence  is  a  process  algorithm 
ordering  that  guarantees  an  accurate  simulation  of  the 
system.  The  ordering  of  process  algorithms  to  form  simula¬ 
tion  sequences  is  the  primary  subject  of  this  dissertation. 

1.2  Literature  Survey 

Most  simulation  languages  contain  an  executive 
routine  that  calculates  simulated  time  and  schedules  appro¬ 
priate  event  routines.  A  literature  survey  of  executive 
routines  of  simulation  languages  has  been  previously  published 
by  this  writer  [1].  This  survey  is  summarized  below. 

For  the  purpose  of  discussion,  simulation  languages 
will  be  divided  into  two  classes,  discrete  event  languages 
and  block  diagram  languages.  A  discrete  event  language  is 
usually  an  extension  of  a  high  level  programming  language. 

For  example,  SIMSCRIPT  is  an  extension  of  FORTRAN.  In  a 
discrete  event  language  a  programmer  defines  the  simulation 
problem  by  specifying  the  state  of  the  simulated  system  and 
subroutines  for  each  change  of  state. 

The  block  diagram  languages  require  no  programming. 
In  these  languages  a  collection  of  basic  building  blocks  is 
supplied  with  the  language.  An  engineer  defines  the  simula¬ 
tion  problem  by  listing  the  blocks  and  the  interconnections 
between  blocks. 

The  problem  of  simultaneous  events  was  initially 
ignored.  It  was  presumed  that  the  order  of  execution  of 
simultaneous  events  would  result  in  little  or  no  difference 
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since  the  simulated  time  step  and  hence  the  resulting  state 
changes  are  usually  small.  In  a  comparison  of  two  discrete 
event  languages,  GPSS  II  and  SIMSCRIPT,  Murphy  [2]  reports 
that  the  two  programs  produced  identical  output  for  a  chosen 
problem  after  adjustments  had  been  made  in  random  number 
generation,  number  rounding,  the  order  in  which  events  are 
caused  and  the  treatment  of  simultaneous  events.  SIMSCRIPT  [3] 
provides  no  aid  in  the  handling  of  simultaneous  events. 

GPSS  [4]  uses  a  chain  propagation  technique,  that  is,  if  an 
event  triggers  a  subsequent  event,  the  subsequent  event  is 
selected  as  the  next  event.  The  SOL  [5]  simulation  language 
employs  a  priority  technique.  The  SIMULA  [6]  simulation 
language  has  a  "quasi-parallel"  operation  mode.  In  this  mode 
a  collection  of  programs  conceptually  operates  in  parallel. 

In  summary,  the  discrete-event  languages  place  the  burden  of 
the  accurate  simulation  of  simultaneous  events  on  the  pro¬ 
grammer  . 

In  the  block  diagram  languages  the  selection  of 
simultaneous  events  is  made  by  the  simulation  program.  It 
will  select  the  order  of  execution  of  behavior  algorithms 
for  the  system  blocks.  The  DAS  [7]  simulation  language 
executes  behavior  algorithms  for  system  blocks  in  the  order 
in  which  block  definitions  are  read  into  the  program. 

Engineers  using  DAS  discovered  that  changing  the  order  of 
block  definitions  resulted  in  different  solutions  for  the 
same  problem.  This  phenomenon  resulted  in  such  a  low  level 
of  confidence  in  DAS  solutions  that  DAS  was  used  primarily 
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to  obtain  scaling  information  for  analog  computers. 

MIDAS  [8],  a  successor  of  DAS,  circumvents  the  above  problem 
by  sorting  the  block  definitions  into  a  predetermined  order. 

The  sorting  of  block  definitions  results  in  an 
ordering  of  simultaneous  events.  The  MIDAS  sort  routine 
rearranges  the  blocks  so  that  all  blocks  have  only  "known" 
inputs  at  all  times.  Inputs  are  "known"  if  they  are  con¬ 
stants,  initial  conditions,  or  outputs  of  previously  ordered 
blocks.  In  MIDAS  each  block  has  a  single  output.  If  all 
outputs  cannot  be  made  "known",  then  an  algebraic  loop  exists 
in  the  block  diagram.  An  algebraic  loop,  also  called  a 
zero  delay  loop,  is  an  instantaneous  cyclic  dependency  of 
block  output  variables  on  block  output  variables. 

An  algebraic  loop  cannot  be  made  "known"  since 
the  cyclic  dependency  requires  that  in  order  to  make  an  out¬ 
put  "known",  it  must  have  been  previously  "known."  MIDAS 
will  accept  an  algebraic  loop  only  if  an  implicit  element 
has  been  specified.  In  such  a  case  the  output  corresponding 
to  the  implicit  element  is  made  "known,"  and  the  sorting 
operation  Is  continued. 

The  BLODI  [9]  language  uses  a  similar  sorting 
technique.  BLODI,  like  MIDAS,  allows  only  one  output  per 
block.  The  BLODI  sorting  program  assigns  one  binary  storage 
cell  for  each  output.  The  two  states  of  each  output  are 
called  "full"  and  "empty."  Initially  all  outputs  which 
represent  inputs  to  delays  are  "full"  and  all  others  are 
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"empty.”  A  block  can  be  compiled  whenever  all  its  inputs 
are  "full"  and  its  output  is  "empty."  When  a  block  is  com¬ 
piled,  its  output  is  marked  "full."  Compilation  continues 
until  all  blocks  have  been  compiled  or  until  no  uncompiled 
block  meets  the  requirements.  In  the  latter  case  the  remark 
"-closed  loop  with  no  delays  or  all  delays"  is  printed  and  the 
problem  abandoned.  Such  a  block  diagram  is  called  a  non- 
admissible  circuit.  The  authors  [9]  state  "there  is  no  way  to 
interpret  a  block  diagram  corresponding  to  a  nonadmissible 
circuit.  To  be  sure  one  could  connect  physical  boxes  in 
such  a  manner  and  something  would  happen.  The  analysis, 
however,  would  involve  the  precise  transient  behavior  of 
the  devices  within  the  pulse  width  -  information  which  is  not 
available  to  the  compiler." 

In  summary,  both  MIDAS  and  BLODI  sort  the  block 
definitions  in  such  a  manner  that  each  block  uses  only  con¬ 
st  anu  or  previously  calculated  values.  The  reader  may  visualize 
this  ordering  as  follows: 

1.  Suppose  that  each  output  has  a  small  fixed  amount 
of  delay  associated  with  it. 

2.  Suppose  that  the  block  diagram  is  quiescent  and 
that  an  input  change  occurs. 

Due  to  the  small  amount  of  delay  associated  with 
each  output,  the  input  change  cannot  propagate  instantaneously 
throughout  the  block  diagram  but  will  ripple  through  the 
block  diagram.  If  the  block  diagram  is  acyclic,  that  is. 
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contains  no  cycles, then  each  block  will  reach  a  quiescent 
state.  The  blocks  that  have  only  block  diagram  inputs  will 
reach  a  quiescent  state  after  a  single  delay  unit.  The  other 
blocks  will  reach  a  quiescent  state  one  delay  unit  later  than 
the  latest  input.  The  order  in  which  the  blocks  reach  a 
quiescent  state  is  essentially  the  same  as  the  order  produced 
by  the  sorting  routine.  It  should  be  noted  that  if  two  or 
more  blocks  reach  the  quiescent  state  at  the  same  time,  the 
order  in  which  these  blocks  are  sorted  is  immaterial.  This 
block  ordering  is  called  the  simulation  sequence. 

If  the  block  diagram  contains  a  cycle, blocks  in 
the  cycle  will  never  reach  a  quiescent  state  and  hence  can 
never  be  ordered. 

Both  MIDAS  and  BLODI  assume  that  each  block  con¬ 
tains  only  a  single  output.  As  a  consequence  of  this  assump¬ 
tion  each  block  will  be  included  exactly  once  in  the  simula¬ 
tion  sequence. 

The  problem  of  constructing  a  simulation  sequence 
for  multiple  output  blocks  has  been  studied  by  Parnas  [10], 
Parnas  defines  a  multiple  output  block,  called  an  augmented 
process,  as  a  5-tuple  (I,  ft,  ip,  Vt  A)  where  E  is  a  set  of 
input  variables,  ft  is  a  set  of  output  variables,  f  is  a  set 
of  state  variables,  y  is  a  single  algorithm  composed  of  a  set 
of  event  description  algorithms,  and  A  is  the  instantaneous 
dependency  relation  of  process  output  values  on  process  input 
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values.  An  interconnection  of  multiple  output  blocks,  called 
a  network,  is  a  quadruple  (A,  Z,  r,  R)  where  A  is  a  set  of 
processes,  Z  is  a  set  of  network  inputs,  ?  is  a  set  of  inter¬ 
connector  variables,  and  R  a  set  of  statements  about  the  net¬ 
work.  R  relates  process  outputs  with  interconneetor  variables. 
Interconnector  variables  with  process  inputs,  and  network 
inputs  with  process  inputs. 

Paruas  considers  two  types  of  algorithm  applications 
called  process  references.  The  first  denoted  by  <A>  is  the 
execution  of  the  process  algorithm.  The  second  denoted  by 
,v<A>  is  the  execution  of  the  process  algorithm  followed  by 
the  restoring  of  the  values  of  ^  to  the  values  they  held 
before  execution.  An  application  of  an  algorithm  will  cal¬ 
culate  a  correct  value  for  an  output  variable  only  if  all 
input  variables  on  which  this  output  depends  have  correct 
values.  A  correct  value  is  obtained  for  the  state  variables 
only  if  all  input  variables  have  correct  values. 

The  simulation  state  of  a  network  is  a  Boolean 
vector  with  one  entry  corresponding  to  each  process  in  A 
and  one  entry  corresponding  to  each  member  in  T.  All  entries 
are  initially  "false."  An  entry  corresponding  to  an  element 
of  r  is  "true"  if  the  interconnector  variable  has  a  correct 
value.  An  entry  corresponding  to  an  element  of  A  is  "true" 
if  all  process  state  variables  have  correct  values.  A  simu¬ 
lation  sequence  is  a  string  of  process  references.  The 
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problem  of  constructing  a  simulation  sequence  for  multiple 
output  blocks  is  that  of  construction  a  simulation  sequence 
with  all  "true"  final  simulation-state  entries. 

Parnas  constructs  a  minimum-cost  simulation  seqr  ice 
from  a  simulation  graph.  The  nodes  of  the  simulation  graph 
are  the  simulation  states  of  the  network.  Each  branch  of  the 
simulation  graph  is  labeled  with  a  process  reference.  A 
branch  depicts  the  change  in  the  simulation  state  that  results 
from  the  application  of  a  process  algorithm.  The  label  of  any 
path  in  this  graph  from  the  ''all-false’'  node  to  the  "all-true” 
node  is  a  correct  simulation  sequence. 

Parnas  proves  in  his  Theorem  3  that  there  exists  an 
"all-true"  node  in  the  simulation  graph  if  and  only  if  there 
exists  no  cyclic  interconnector  dependency.  In  other  words, 
a  correct  simulation  sequence  exists  if  and  only  if  the  net¬ 
work  contains  no  zero-delay  loop 

In  his  conclusion  Parnas  states:  "It  appears  quite 
easy  to  generate  networks  in  which  the  graph  used  in  searching 
for  optimal  solutions  becomes  too  large  to  handle." 

1.3  Approach  to  Solution  of  Problem 

The  formulation  of  the  problem  is  essentially  similar 
to  that  of  Parnas  [10].  Functional  dependencies  will  be  repre¬ 
sented  by  relations.  Only  the  execution  of  a  given  set  of 
algorithms  will  be  considered.  Equivalent  sequences  involving 
the  application  of  other  algorithms  will  not  be  considered. 
Interconnector  variables  will  not  be  employed.  The  single- 


algorithm  per  process  assumption  is  replaced  by  a  set  of 
algorithms  per  process  assumption.  A  single  algorithm  in 
each  process- algorithm  set  is  assumed  to  update  the  set  of 
process  state-variables.  This  assumption  will  lead  to  a 
normal  form  for  simulation  sequences  and  will  simplify  the 
-formation  of  simulation  sequences  by  removing  the  dependency 
of  simulation  sequences  on  state  variables.  It  is  shown  in 
section  8.1  that  this  assumption  produces  no  loss  of  generality. 

The  process  input-output  dependency  relation  is 
composed  with  the  network  interconnection  relation  to  form  a 
network  output- dependency  relation.  This  output-dependency 
relation  is  represented  by  a  directed  graph.  The  graph  is 
acyclic  since  no  zero-delay  loops  are  allowed  in  the  network. 

The  theory  of  directed  graphs  will  be  applied  to 
obtain  a  minimum- length  output  simulation  sequence.  The 
strategy  employed  is  partially  to  order  the  set  of  algorithms. 

An  algorithm-dependency  relation  is  formed  by  condensing 
the  transitive  closure  of  the  output-dependency  relation  ac¬ 
cording  to  the  algorithm  definition.  This  algorithm-de¬ 
pendency  relation  need  not  be  acyclic. 

If  the  algorithm-dependency  relation  contains 
reflexive  elements,  that  is,  if  there  exists  an  algorithm  with 
two  or  rrore  related  outputs,  this  algorithm  is  split.  A  new 
algorithm-dependency  relation  is  formed  using  the  modified 
algorithm  definition.  Although  this  relation  is  irreflexive 
it  need  not  be  acyclic.  An  ordering  of  the  modified  algorithm 
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set  is  generally  not  possible.  An  ordering  of  the  maximal 
connected  subsets  is,  however,  always  possible.  Furthermore 
a  concatenation  of  minimum-length  simulation  sequences  for 
maximal  connected  subsets  according  to  their  ordering  is 
shown  to  form  a  minimum-length  simulation  sequence. 

The  formulation  of  the  problem  is  discussed  in 
Chapter  II.  In  Chapter  III  the  properties  of  simulation 
sequences  are  derived.  A  normal  form  is  developed  and 
the  problem  of  constructing  a  simulation  sequence  is 
reduced  to  the  problem  of  constructing  a  simulation  sequence 
for  output  variables.  It  is  shown  that  state-variable 
dependencies  can  be  ignored.  Level  partitions  of  the  set  of 
output  variables  are  investigated  in  Chapter  IV.  Algorithm 
ordering  is  investigated  in  Chapter  V.  Algorithm  splicting 
is  investigated  in  Chapter  VI.  The  formation  of  simulation 
sequences  is  discussed  in  Chapter  VII.  Extensions  to  a 
larger  class  of  problems  are  discussed  in  Chapter  VIII. 
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II.  SYSTEM  MODEL 


The  system  to  be  simulated  is  assumed  to  be  char¬ 
acterized  as  a  network  of  processes. 

2 . 1  Process  Definition 

Each  process  must  be  capable  of  modelling  a  sub¬ 
system.  A  basic  system  model  is  shown  in  Figure  1.  In  this 
model  a  system  consists  of  a  set  of  input  variables,  Z;  a 
set  of  output  variables,  Y:  a  set  of  state  variables,  M;  and 
a  set  of  memory-free  functions,  H.  For  the  purpose  of  simula¬ 
tion,  this  model  is  deficient  in  the  following  respects: 

1.  The  model  implies  that  the  value  of  each  output 
or  state  variable  depends  upon  the  value  of  every 
input  and  state  variable.  Although  this  assumption 
causes  no  theoretical  difficulty,  it  places  a 
severe  restriction  on  the  problem  of  obtaining  an 
accurate  simulation  sequence  for  the  system.  In 
general,  an  output  or  state -variable  value  is  only 
dependent  upon  some  input  and  state -variable  values. 
The  knowledge  of  the  functional  dependency  of  process- 
variable  values  will  allow  the  formation  of  an  ef¬ 


ficient  simulation  sequence.  The  functional  depend¬ 
ency  is  specified  by  the  R  relation.  Process  variables 
i  and  j  are  R- related,  denoted  by  (i,j)  C  R  or  iRj  , 
if  the  value  of  variable  j  depends  upon  the  value  of 


variable  i. 


2.  For  the  purpose  of  simulation,  it  is  often  desirable  | 

to  combine  several  functions  into  a  single  algorithm.  | 

When  functions  are  combined  into  an  algorithm,  it  | 

is  frequently  possible  to  use  their  similarities  to  | 

I 

reduce  the  amount  of  computation  required  by  elimin-  | 

ating  duplicate  calculations.  The  combination  of  | 

I 

functions  into  algorithms  is  specified  by  the  A  | 

I 

(algorithm)  partition.  Initially  we  will  assume  1 

that  A  partitions  YUM.  Later,  in  order  to  guarantee  I 

that  an  accurate  simulation  exists  for  every  process, 
the  A  partition  will  be  further  restricted. 

Formally, a  process  p  is  defined  to  be  a  quintuple 

p  =  (2),Y,M,R,A) 
l  -  finite  set  of  inputs 
Y  -  finite  set  of  outputs 
M  -  finite  set  of  state  variables 
R  C  (I  U  M)x(Y  UM) 

A  -  partition  of  (Y M) 

The  sets  of  input,  output,  and  state  variables  are 
assumed  to  be  disjoint,  i.e., 

ZflY  =  I  flM  =  YOU  =  $  (empty  set) 

Capital  letters  will  be  used  to  denote  sets,  and 
small  letter  will  be  used  to  denote  elements  of  sets.  In 
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circumstances, where  there  exists  a  need  to  differentiate 
between  members  of  a  set,  subscripts  will  be  used  to  denote 
specific  elements.  Subscripts  will  be  used  on  sets  to  denote 
subsets.  Thus  a  £  Z  and  Z^C  Z.  Superscripts  will  be  used 
to  denote  the  process  to  which  an  element  or  a  set  belongs. 

Each  input,  output  and  state  variable  has  a  value 
as  an  attribute,  l.e., 


v^  -  value  of  variable  i 

An  ordered  set  of  variables,  Z,  has  a  value  vector, 
Vz,  as  an  attribute,  i.e.. 


-  v.  ,v.  , . . . , v. 

X1  12  V 


where 


v.  -  value  of  variable  i .  £  Z 
j  J 


Corresponding  to  each  output  and  state  variable 
there  exists  a  function,  h,  which  determines  its  variable  value, i.e.. 


i  =  hi(VZi,VM  r°r  each  1  6  Y  U 


where  Zi  =  {a|crRi} 
=  {m|mRi} 


and 


.  u 


The  R  relation  determines  the  domain  of  each 
function.  For  the  purpose  of  discussion,  it  will  be  con¬ 
venient  to  differentiate  between  output  functions  and  state- 
variable  functions. 

f  -  output  function 
g  -  state-variable  function 

The  set  of  output  functions  and  the  set  of  state 
variable  functions  are  denoted  by  F  and  G. 

F  =  {f  }  =  ih±\l  e  Y} 

G  =  (g  }  =  (h1|i  €  M} 

By  definition  there  exist  one-to-one  mappings 

between  F  and  Y  and  between  G  and  M.  Hence  h^  will  often  be 
denoted  by  i.  This  partition  of  the  set  of  functions  also 
partitions  the  domain  of  the  functions  and  hence  specifies 
a  partition  of  R. 

R  =  RAu  rburcurd 
R-a  C  MxM 
rb  C  zxm 
Rc  C  MxY 


SP£K«?!S3=?  055*j*#k^«=- 
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These  process  definitions  are  displayed  in 

Figure  2. 

In  a  simulation  problem,  separate  algorithms  may 
not  exist  for  each  function.  For  various  reasons,  several 
functions  may  have  been  lumped  into  a  single  algorithm.  The 
A  partition  specifies  which  functions  are  to  be  executed 
concurrently.  Although  A  is  specified  as  a  partition  of 
(YUM),  it  will  be  interpreted  as  a  partition  of  (FUG). 

A  =  {H^} 

H1  n  =  $  for  i  ft  j 

U  H  =  p  U  g 
i  1 

As  a  consequence  of  the  A  partition,  each  output  and 
state-variable  function  is  an  element  of  exactly  one  algorithm. 
All  algorithms  are  nonempty  sets  of  functions.  An  empty  set 

of  functions  will  not  change  any  variable  values  and  hence  is 
not  of  interest. 

An  algorithm  sequence  a  is  an  ordered  string  of 
algorithms.  X  is  the  algorithm  sequence  of  length  zero. 
Superscripts  are  used  on  value  vectors  to  denote  the  appli¬ 
cation  of  an  algorithm  sequence. 


i.ttw  ’3gS 


FIGURE  2 
PROCESS  MODEL 


i «  =_ 

4.  V 


-  initial  value  vector 

Va  -  value  vector  after  the  application  of 
algorithm  sequence  a. 

In  the  application  of  an  algorithm  all  functions 
in  the  algorithm  are  assumed  to  be  applied  in  parallel,  i.e., 
all  functions  use  the  value  of  variables  prior  to  the  appli¬ 
cation  of  the  algorithm.  Thus  if  the  algorithm  H  is  applied 
after  the  algorithm  sequence  a,  the  resulting  values  are: 


where  i  €  (Y  U  M)  ^v“  otherwise 

The  reader  will  note  that  the  algorithm  H  only  alters  the 
values  of  variables  corresponding  to  functions  in  H.  All 
functions  in  H  use  the  values  of  variables  prior  to  the 
application  of  H.  This  differs  from  a  sequential  application, 
in  which  the  initial  values  are  not  saved  and  the  values 
resulting  from  the  functional  calculation  are  employed  in 
subsequent  calculations. 

A  simulation  sequence  for  a  process  is  any 
algorithm  sequence  which  calculates  new  values  for  each 
output  and  state  variable  using  only  the  initial  values  of 
the  input  and  state  variables.  A  simulation  sequence  is 
formally  defined  in  terms  of  the  simulation  state. 


-  19  - 


The  simulation  state  of  a  variable  is  an  attribute 
determined  by  the  variable  value.  Initially  the  simulation 
state  of  each  output  and  state  variable  is  equal  to  "0".  The 
simulation  state  of  a  variable  is  equal  to  ”1”  if  a  correct 
value  has  been  calculated  for  this  variable,  and  the  simulation 
state  is  equal  to  "2"  if  the  variable  value  is  incorrect. 

si  -  simulation  state  of  variable  i  €  Zlj  Y  (JH 


s 


i 


0 


where 


1 


if  YUM 


l2 


if 

if 


v_. 


otherwise 


The  simulation  state  of  a  process  Input  variable  is 
determined  by  the  network  structure.  It  is  equal  to  the 
simulation  state  of  the  network  variable  determining  its 
value . 

The  simulation  state  of  an  ordered  set  of  vari¬ 
ables  is  represented  by  the  simulation-state  vector,  denoted 
by  a  capital  S.  This  vector  as  the  value  vector  is  a  func¬ 
tion  of  previously  applied  algorithms. 


S»SS WST¥ 
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S^  -  simulation  state  vector  with  a  single  component 
associated  with  each  element  of  Z 
-  initial  simulation  state  vector 
Sa  -  simulation  state  vector  after  the  application 
of  algorithm  sequence  a. 

In  the  above  definition,  the  simulation  state 
depends  upon  the  variable  value  and  hence  upon  the  variable 
function.  It  is  conceivable  that  a  correct  value  may  be 
calculated  for  a  variable  although  incorrect  values  are  used 
for  the  variables  on  which  this  variable  depends.  The  above 
situation  is  removed  by  defining  the  simulation  state  in 
terms  of  a  initial  simulation  state  and  rules  which  do  not 
explicitly  depend  upon  variable  values. 

The  application  of  an  algorithm  changes  the  simu¬ 
lation  state  via  the  following  rule. 

if  h.  H 

if  (h^  C  H)  and  aRi  =>  s0  =  lj 
and 

otherwise 

An  algorithm  H  will  not  alter  the  simulation  state  of  variables 
which  are  not  contained  in  H.  If  the  variable  is  included  in 
H,  the  simulation  state  is  set  equal  to  "l"  iff  all  R-related 


mRi  =>  s_ 


=  0 

4 
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input  variables  have  "1"  simulation  state  and  all  R  related 
state  variables  have  "0"  simulation  state.  If  the  variable 
is  included  in  H  and  the  above  conditions  are  not  satisfied, 
the  simulation  state  is  set  equal  to  "2". 

A  simulation  sequence  for  a  process  is  any  algorithm 
sequence  a  such  that 


si  =  0  for  each  i  €  Y  (J  M 


si  =  l  for  each  i  e  £ 


s T  =  1  for  each  i  €  JU  YUM 


A  simulation  sequence  does  not  exist  for  all  pro¬ 
cesses.  Consider,  for  example,  the  process  shown  In  Figure  3< 
For  this  process 

H1  =  {m1,y} 


H2  =  {m2) . 


No  simulation  sequence  exists  for  this  process  since 

HiH2 

s  1  =  2 

m2 


H2H1 


2. 


_  oo 

—  k  W  — 


P  «  ({o},{y},  {mltm2),  R, 

R:  (a,  n^),  (in^nig),  (m2,y),  (m2,m1) 
# 

Hi  “  {mi»y> 

H2  * 


FIGURE  3 

Process  For  Which  No  Simulation  Sequence  Exists 


If  H1  is  executed  first,  a  new  value  is  computed 
for  m1.  H2  will  use  the  new  value  for  and  hence  will 
calculate  an  incorrect  value  for  m2.  If  H2  is  executed  first, 
the  value  of  calculated  by  H±  is  incorrect.  One  solution 
to  this  dilemma  is  to  store  the  initial  value  of  all  state 
variables,  and  only  use  these  stored  values  in  all  calcula¬ 
tions.  This  technique  requires  double  storage  for  all  state- 
variables.  An  alternative  solution  is  to  lump  all  state 
variable  functions  into  a  single  algorithm.  In  this  thesis, 

it  will  be  assumed  that  the  A  partition  generates  one  algorithm 
equal  to  G,  i.e.. 


A  -  {F±}  U  G 

FinV  * 

U  F.  =  P 
i  1 


Under  this  assumption  a  simulation  sequence 
exists  for  every  process.  The  reader  will  note  that 


a  =  F  F0F0. .  .f  g 
12  3  n 


where 


n 

U 

k=l 


F 


} 

t 

l 


is  a  simulation  sequence  for 


any  process. 


X&ft&tor VJA&ij'.Vr  •  •*  .{.*»,  i>.t  ays*  *v.*  ■ 
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The  simulation  state  and  hence  a  simulation  sequence 
is  independent  of  variable  functions  and  variable  values.  Due 
to  the  one-to-one  correspondence  between  functions  and  output 
and  state  variables;  an  algorithm  may  be  represented  by  a  set 
of  output  and  state  variables.  This  representation  of  algorithms 
will  be  used  in  the  remainder  of  this  thesis.  The  symbols  P,  G, 
and  H  will  still  be  used  to  denote  algorithms  to  avoid  confusion 
between  algorithms  and  variable  sets. 

2.2  Network  Definition 

A  network  is  a  collection  of  processes.  In  addition 
to  a  list  of  processes,  the  network  definition  must  specify  the 
interconnections  between  processes.  This  is  accomplished  by 
the  RN  relation.  It  relates  the  entire  set  of  process  outputs, 
denoted  by  Yp,  to  the  entire  set  of  process  inputs,  denoted  by 
Ip.  The  relation  also  specifies  the  connection  of  network 
inputs  to  process  inputs.  The  network  outputs  Y^  specify 
the  simulation  output.  The  inclusion  of  network  outputs  in 
the  network  definition  allows  the  network  definition  to  serve 
as  a  process  definition.  Formally,  a  network  N  is  defined  to 
be  a  quadruple 

P  =  {p}  -  finite  set  of  processes 
IN  -  set  of  network  inputs 
Y>j  9  Yp  -  set  of  network  outputs 
R^  9  (Yp  u  Z^)  *  Zp  -  network  relation. 
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A  network  is  thus  an  interconnection  of  one  or  more 
processes.  Processes  are  assumed  to  have  no  variables  in 
common,  for  example, 

pi  pi 

£  n  Z  J  =  $  (empty  set)  for  i  /  j  , 

The  set  consisting  of  the  union  of  all  process  sets  of  a 
given  type  will  be  denoted  by  the  subscript  P. 

£p  =  U  £p 
P 

Network  inputs  represent  external  forces  on  the 
system.  The  value  of  a  network  input  variable  is  assumed 
to  be  given  for  all  time.  A  network  output  is  a  directly 
observable  process  output. 

The  network  relation  R^  specifies  the  connection  oi 
process  outputs  and  network  inputs  to  process  inputs. 

(y,a)  €  Rn  or  yRNa  if  process  output  y  is  con¬ 
nected  to  process  input  a. 

(aN,a)  c  Rn  or  if  network  input  aN  is 

connected  to  process  input  a. 

A  network  is  complete  if  every  process  input  is 
connected  to  a  single  process  output  or  network  input. 


‘  f  -w-i 


1 
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Formally  a  network  is  complete  if  for  each  o  €  Ep  there  exists 
a  unique  i  such  that  iR^a.  The  reader  will  note  that  the 
behavior  of  an  incomplete  network  is  not  defined.  In  the 
following  a  network  is  assumed  to  be  complete. 

The  network  relation  specifies,  for  each  process 
input,  the  process  output  or  network  input  which  determines 
its  value.  The  simulation  state  of  a  process  input  is 
equal  to  the  simulation  state  of  the  process  output  or  net¬ 
work  input  to  which  it  is  connected,  that  is,  for  each 


(i»cr)  €  sa  =  s^. 


A  simulation  sequence  for  a  network  is  any  algorithm  sequence 
a  such  that 

s^  =  0  for  each  i  e  Yp  [J  Mp 


and 


s* 

i 


=  1 


for  each 


1  €  Z 


N 


s“  =  1  for  each  i  €  Yp  UMp. 
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2 . 3  Example  Problem 

Simulation  sequences  do  not  exist  for  all  network 
Before  discussing  the  properties  of  a  simulation  sequence, 
a  network  example  is  presented.  A  simulation  sequence  will 
be  constructed  for  the  network  of  Figure  ^4. 


R  ( { Pj >?2 * P g ^ 1 ^ ^ 1  * y 7 ^ 1 ^ 

"  ^°N ,a  2  ^  >^y2,05^  »  (y  ^  *^y4*^4^  * 

(y4,tf7),(y5,a8),(y6,a3) 

P-^  -  3 ^  *  ^1  3 *  ^  *  ^^1  ^  ^ ^ ^1  *^2  3 ^  *  ^^1^ 

R:  (a1,y3),(a1,m1),(a2,m1),(a3,y1),(m1,y1),(m1,y2) 

P2  =  1 1 0]j  ,Cg } ,  {y^  ,y ^ } ,  {m2 ,m3 }  ,R,^{y|j  ,y^ } , {m2  ,m3}|’ j 

R:  (0]j,y5)  »(a5,y5)5  (o5,m2),(a6,m2),(m2,m3),(m2,y2l), 
(m3  >y^ ) > (m3  >y 5 ) 

P3  =  {a7,a8,a9},{y6,y7},{m1|},R,|{y6,y7},{m2j}}j 

R •  (o7,yg),(og,y7),( Pq  >  y  7 ) > ( » ^4  ^  ^  ^4  *  ^4  ^  ^  ^4  * y  7  ^ 


The  above  definitions  are  displayed  in  Figure  4. 

The  R^  relation  is  represented  by  solid  lines  and  the  R 
relations  are  represented  by  broken  lines.  The  R  relation  is 
defined  to  be  the  union  of  the  process  relations.  The  R 


ft-  -” 
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relation  is  thus  between  (IpU  Mp)  and  (MpU  Yp).  The 
composition  of  R^  and  R,  written  R^R,  is  defined  as  follows 

w(R^R)z  <=>  there  exists  a  €  Ip  such  that 
wR^o  and  aRz 

The  R^R  relation  is  between  (1^  UYP)  and  (MpU  Yp).  The 
RNR  relation  specifies  the  dependency  of  process  output  and 
state  variables  on  network  input  and  process  output  variables. 
For  the  above  example 

RjjF*  ) » (y2  ^  *  ^2  *^5  ^  * 

(y  3 ) * (y  3 *^4  ^,^Y3*y7^»^Y^»yy)> 

(y4  ,y5 ) ,  (y;, ,!«! ) » (y  4  »y  3) ,  (y^  >y6 )  > 

(yg>y1)* 

The  predecessor  of  an  element,  z,  in  the  range  of 
a  relation  R,  written  R(z),  is  defined  to  be 

R(z)  =  {w|wRz}. 

In  the  network  example 

R(y5)  =  {o4,a5,m3} 

RjjR(y^)  -  {y2»y4>- 


and 
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The  eligible  set  E  is  the  set  of  process  output 
and  state  variables,  whose  simulation  state  is  not  equal 
to  1,  but  whose  simulation  state  can  be  set  equal  to  1  by 
the  application  of  a  single  process  algorithm. 

E  =  jv](sv  5 *  1)  and  (sv/  =  1  for  the  v#  €  R^R(v)) 

and  (sm  =  0  for  the  m  £  R(v))j 

The  eligible  set  is  a  function  of  the  simulation  state.  For 
example,  the  simulation  state  vector  S  is  defined  as  follows: 


The  simulation  state  of  aN  is  always  equal  to  1  and  is  not 
included  in  S.  The  intial  simulation  state  vector  is 

SX  =  (0000000  0000) 

and  for  this  simulation  state 

E  *  •fy2>y4^3>- 

The  initial  choice  of  process  algorithms  is  thus  reduced  to 

p 

three  choices.  If  G  is  selected,  it  will  be  impossible  to 
calculate  correct  values  for  y^  and  y,-  since  the  initial  value 

p 

of  m2  and  m^  may  be  altered  by  the  G  =  {m^m^}  algorithm. 

A  safe  criteria  for  the  selection  of  process  algorithms  is 
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to  update  state  variables  last.  If  a  minimum -length  Simula- 

p 

tion  sequence  is  desired,  F  should  be  selected.  The  reason 
for  selecting  F^  over  F^  is  not  obvious  and  will  be  developed 

p 

in  succeeding  chapters.  Selecting  (y^y^)  s  F  as  the  initial 
algorithm  will  calculate  a  correct  value  for  y^  and  an  in¬ 
correct  value  for  y^. 

•ci2 

S  =  (0001200  0000) 

The  application  of  F^  removes  y^  from  E  and  adds  y^,  yg  and 
m1  to  E. 

E  =  {y  2  »y  2  *^i  ^ 

l  o  i  2 

The  choice  of  algorithms  is  thus  between  F  ,  F  ,  G  and  G  . 

*i  p 

The  selection  of  G  or  G  will  deny  the  formation  of  a 
simulation  sequence.  F1  will  correctly  update  two  outputs, 
and  F-5  will  correctly  update  one  output.  Intuitively,  F 
should  thus  be  preferred  over  F^  if  a  minimum— length  simula¬ 
tion  sequence  is  desired.  Actually,  the  number  of  outputs 
correctly  updated  is  a  poor  criteria  for  the  selection  of 

process  algorithms.  In  this  case,  either  selection  can  lead 

3 

to  a  minimum- length  simulation  sequence.  We  will  select  F  . 

2  3 

SF  F  =  (0001212  0000) 

and 


E  =  {y1,y2,y3,m1,m3} . 
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We  are  forced  to  select  F1 


p2p3pl 

S  =  (1111212  0000) 


^  —  <y5>  ^  i®  2  *  ^  * 

All  the  outputs  of  p^  have  1  simulation  state 
entries.  G1  can  thus  be  selected.  G2  and  G3  cannot  as  yet 
be  selected.  Our  choices  are  thus  F2  and  G1.  Either  choice 
can  lead  to  a  minimum-length  sequence.  We  select  G1. 


p2„3_l  l 

S  =  (1111212  1000) 


1-1  j^2  >^2  J  • 


We  are  forced  to  select  F2. 


W  2  p  3  p  1  n  1 TT1  2 

SP  =  (1111112  1000) 


G3  should  still  not  be  selected.  F3  and  G2 

o 

We  select  F. 


are  equal  choices. 


and 
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„2„3plrlp2„3 

SF  F  F  G  P  F  =  (limil  1000) 


E  =  {m2  jin^  ,m^} 


G2  and  G3  are  equal  choices.  We  select  G3. 


„2„3„lpl„2„3p3 

SF  F  F  G  F  F  G  =  (1111111  100 


and 


E  =  {m^m^} 


Only  G  remains 


sFWgVpW  _  anlnl  lm) 


F2F3F1G1F2F3G3G2  is  a  minimum-length  simulation  sequence, 
Other  minimum-length  simulation  sequences  are 


F2F3plF2F3G2G3Gl 


F2F^F1G1F2G2F^G^ 


F2F1F2G2F3G3F1G1 


F2FlF2F3p1GiG2GJ< 
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If  F^  is  selected  as  the  initial  algorithm,  simu¬ 
lation  sequences  are 


f1f2f2f1f3g3f1g2g1 


F1F2F3F1G1F3G3F2G2. 


Although  many  simulation  sequences  exist  for  the 
example  network,  there  are  many  networks  for  which  no 
simulation  sequence  exists.  For  example,  if  (cJgjy^)  is 
added  to  the  R  relation  in  the  example  problem,  no  simulation 
sequence  will  exist  for  the  network.  The  addition  of  the 
(c^y^)  term  creates  a  cyclic  chain  of  input-output  and  output- 
input  dependencies,  namely 

(a6,yi,)(yi|,a1)(a1,y3)(y3,a6). 

Due  to  the  cyclic  nature  of  this  dependency  chain,  in  order 
for  a  variable  in  the  chain  to  achieve  ”1"  simulation  state 
it  must  have  previously  had  "1”  simulation  state.  Since  the 
initial  simulation  state  of  an  output  variable  is  always 
equal  to  "0",  no  output  in  the  chain  can  achieve  a  "1" 
simulation  state.  Hence  no  simulation  sequence  will  exist. 
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III.  SIMULATION  SEQUENCE  PROPERTIES 

Although  it  is  possible  to  develop  a  large  number 
of  simulation-sequence  properties,  only  those  properties  which 
provide  a  clue  of  the  length  or  structure  of  a  minimum-length 
simulation  sequence  are  developed  in  this  chapter.  In  Section 
3.1  properties  which  follow  as  a  consequence  of  the  simulation 
sequence  definition  are  developed.  Properties  which  depend 
upon  the  structure  of  the  network  are  developed  in  Section  3.2. 

3.1  Basic  Properties 

For  any  simulation  sequence  the  initial  simulation- 
state  vector  always  contains  all  ”0"  entries.  The  final  simulation 
state  vector  contains  all  "1”  entries.  A  simulation-state  entry 
can  only  change  by  the  application  of  an  algorithm.  Since  the 
set  of  algorithms  partitions  the  set  of  network  variables,  each 
algorithm  must  be  applied  at  least  once. 

Property  1  -  Covering 

In  any  simulation  sequence  each  algorithm  is  included 
at  least  once. 

Property  2  -  Minimum  Length 

The  length  of  any  simulation  sequence  is  greater  than 
or  equal  to  the  rank  of  the  algorithm  partition. 

Property  2  places  a  lower  bound  on  the  length  of  a 
simulation  sequence.  An  upper  bound  on  the  length  of  a 
minimum-length  simulation  sequence  will  be  developed  in 
Chapter  IV. 
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Property  3  -  Normal  Form 

Given  any  simulation  sequence  a,  let  cty  denote  the 
algorithm  sequence  obtained  from  a  by  deleting  all  state- 
variable  algorithms,  and  let  aM  denote  an  arbitrary  ordering 
of  the  set  of  state-variable  algorithms.  The  algorithm  se¬ 
quence  otyaM  formed  by  concatenating  cty  and  is  a  simulation 
sequence. 

Proof 

In  order  to  set  the  simulation  state  of  an  output 
equal  to  "l",  the  simulation  state  of  every  state  variable  to 
which  this  output  is  R  related  must  be  equal  to  "0",  and  the 
simulation  state  of  every  output  variable  on  which  this  output 
depends  must  be  equal  to  "l”.  Delaying  the  state-variable 
algorithms  will  guarantee  that  the  simulation  state  of  every 
state  variable  is  equal  to  "0"  during  the  application  of  the 
output  algorithms.  Furthermore  since  the  order  of  execution 
of  output  algorithms  has  not  been  altered,  every  output  variable 
will  have  "1"  simulation  state  after  Oy. 

1o  set  the  simulation  state  of  a  state  variable  to 
"1",  all  process  inputs  to  which  the  state  variable  is  R-related 
must  have  "1”  simulation  state  and  all  state  vs.riables  to  which 
it  is  R-related  must  have  ”0"  simulation  state.  After  the  out¬ 
put  algorithm  sequence  cty  has  been  applied,  all  outputs  have  "1" 
simulation  state.  Network  inputs  always  have  "1"  simulation  state. 
Thus  all  process  inputs  have  "l”  simulation  state  after  cty.  A 
state  variable  can  only  be  R  related  to  a  state  variable  in  the 
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same  process.  The  assumption  of  one  state  variable  algorithm 
per  process  guarantees  that  every  process  state  variable  will 
have  "0"  simulation  state  prior  to  the  application  of  the  state- 
variable  algorithm  corresponding  to  this  process.  The  simu¬ 
lation  state  of  every  output  and  state  variable  is  thus 
equal  tc  "l”  after  ayaM‘  Q.E.D. 

As  a  consequence  of  Property  3,  each  state  variable 
algorithm  is  executed  exactly  once  in  and  the  order  of 
execution  of  state-variable  algorithms  within  is  immaterial. 
The  problem  of  constructing  a  simulation  sequence  is  thus 
reduced  to  the  problem  of  constructing  cty  for  a  simulation  se¬ 
quence.  The  algorithm  sequence  ay  will  be  called  an  output 
simulation  sequence. 

In  Section  2.3  a  simulation  sequence  was  formed  for 
an  example  network  via  a  prudent  choice  of  algorithms.  The 
selection  of  particular  algorithm  sequences  was  shown  to  deny 
the  formation  of  a  simulation  sequence.  This  difficulty  does 
not  occur  in  the  formation  of  output  simulation  sequences.  In 
fact,  if  a  simulation  sequence  is  known  to  exist  for  a  network, 
any  output  algorithm  sequence  can  be  selected  as  the  initial 
segment  of  an  output  simulation  sequence. 

Property  ^  -  "1"  Simulation  State 

If  the  simulation  state  of  an  output  variable  is 
equal  to  "1"  after  some  initial  segment  of  a  finite-rlength 
output  algorithm  sequence,the  application  of  successive  output 
algorithms  will  not  alter  the  simulation  state. 
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Proof  -  By  Contradiction 


Consider  an  output  y  for  which  s  =  1  after  some 

•7 

Qt  P 

output  algorithm  sequence  a.  Assume  y  c  P  and  that  s  /  1. 

v 

Clearly  s  cannot  be  equal  to  ”0”  and  must  hence  be  equal  to 

«7 

"2" .  Thus  there  exists  some  y,  <i  RNR(y)  such  that  sa  /l. 

Since  s  was  previously  equal  to  "l," sa  must  be  equal  to 
yl  yl 

,,2".  Hence  we  may  conclude  that  a  is  of  the  following  form: 


a  =  where  y^  € 


and 


2. 


1 

If  the  length  of  a,  equals  zero  s  =  1  contradicts  the  initial 

X  yl 

condition.  If  is  greater  than  zero,  repeat  the  above 

*  / 
argument  using  y^  instead  of  y  to  obtain  a^.  Since 

/ 

#a ^  this  procedure  must  terminate  and  thereby  contradict 

the  original  assumption.  Q.E.D. 

Property  ^  states  that  once  the  simulation  state  of 
an  output  is  equal  to  "1",  the  application  of  output  algorithms 
will  not  alter  the  simulation  state  of  this  output.  In  the 
construction  of  an  output  simulation  sequence  it  is  thus  only 
necessary  to  achieve  a  "l"  simulation  state  for  each  output  at 
some  point  in  the  algorithm  sequence.  Only  techniques  for 
attaining  a  "l”  output  simulation  state  need  be  investigated. 

3 . 2  Structural  Properties 

In  Section  3*1,  the  basic  properties  of  simulation 
sequences  were  presented.  These  properties  are  a  consequence 
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of  the  definition  of  a  simulation  sequence,  and  are  independent 
of  the  structure  of  the  network  and  process  relations.  In 
this  section  a  criterion  for  the  existence  of  simulation  se¬ 
quences  is  developed.  A  simulation  sequence  exists  for  a 
network  if  ay  exists,  can  always  be  constructed. 

In  the  formation  of  an  output  simulation  sequence, 
only  the  dependency  among  output  variable  values  need  be 
considered.  The  output  dependency  relation  Ry  is  defined 
as  follows : 

ylRYy2  lff  ylRNRy2 
For  the  example  problem  of  Section  2.3 

i 

Ry :  (y2,y5)Ay3,y7)Ay5>y7)  I 

! 

% 

(yn>y5),(yii>y3)>(yii>y6)>(y6>yi')  1 

5- 

5 

The  Ry  relation  is  the  direct  output  dependency 
relation.  In  the  example  the  value  of  yy  depends  directly 
upon  the  values  of  yr  and  y3  since  y^  Ryy7  and  y3RYy7*  rhe 
value  of  y5  depends  upon  the  value  of  y2  and  the  value  of 
depends  upon  the  value  of  y^.  The  value  of  y^  thus 
depends  indirectly  upon  the  values  of  y2  and  y^.  This 
indirect  dependency  relation  is  the  transitive  closure  of  Ry. 

Given  any  relation  R  between  a  set  Y  and  itself,  the 
transitive  closure  of  R,  written  is  defined  as  follows: 
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R+  *  U  R1  =  R  U  (RR)  U  (RRR)  U  ... 
i=l 

If  Y  is  finite  and  contains  m  elements,  only  the 
first  m  terms  need  be  included  in  the  union. 

For  the  example  problem 

Ry  •  (y2»y5My2jy7My3>y7), 

(y4»y5)»(y4»y3)*(y4*yg), 

(y4>y7),(y4,y1),(y6,y1) 

A  relation  R  is  called  acyclic  if  R1"  is  irreflexive, 
i.e.,  R  is  acyclic  if 

for  each  y  €  Y  (y,y)  R+ 

The  relation  R  is  called  cyclic  if  there  exists  some  y  such 
that  yR+y.  The  example  relation  is  acyclic. 

Any  relation  R  defined  on  a  set  Y  has  a  convenient 
graphical  interpretation.  The  elements  of  Y  are  the  nodes 
of  the  graph,  and  the  elements  of  R  are  the  edges  of  the 
graph . 

A  directed  graph  is  a  system  Y  =  <Y,R> 

Y  -  set  of  nodes 
R  -  incidence  relation 
R  C  Yxy 
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The  graph  contains  an  edge  from  y^  to  y2  if 

(yx,y2)  €  R. 

The  network  output  graph  Yp  is  defined  as  follows: 

Yp  *  <Yp,Ry> 

The  network  output  graph  for  the  example  problem  is 
shown  in  Figure  5- 

A  path  from  yQ  to  yn  in  a  graph  Y,  written  P(yQ,yn), 
is  a  sequence  of  edges 

P(y0,yn)  =  ((y0,y3)(yl,y2)*'*(yn-l,yn)) 

such  that 


and 


(yi_l,yi)  €  R  for  all  0  <  i  <  n 


except  for  yfi  and  y^ 


i  /  j  implies  that  yi  /  y^  . 

The  length  of  a  path,  written  #P,  is  the  number  of 
edges  in  the  path.  The  P(y,y)  consisting  of  the  empty  se¬ 
quence  of  edges  is  a  path  of  length  zero. 

A  cycle  is  a  path  of  nonzero  length  for  which  the 
initial  node  and  final  node  coincide.  A  graph  will  contain 
a  cycle  if  R  is  cyclic. 


FIGURE  5 


NETWORK  OUTPUT  GRAPH  FOR  NETWORK 
DISPLAYED  IN  FIGURE  ¥ 


I 

The  output  dependency  relation,  Ry,  is  a  precedence 
relation  concerning  the  "1"  simulat'  r'n  state.  The  reader  will 
note  that  if  outPut  y^  must  have  a  "l”  simulation 

state  before  the  simulation  state  of  y2  can  be  set  equal  to 
"1".  Furthermore  if  yp  €  Fp  and  y2  C  F2>  the  output  algo¬ 
rithm  F.  must  precede  the  algorithm  F2  in  any  simulation  sequence. 
Property  5  -  Path  Ordering 

Given  a  network  output  graph  Yp  =  <Yp,Ry>  and  a  path 

p(y1,yn)  =  ((y1,y2)(y2>y3)>-**»(yn_1>yn)) 

of  Yp,  assuming  yp  c  F^  for  all  1  £  i  ^  n  then  any  output  sim¬ 
ulation  is  of  the  form 


anF,a,F~a0F0. .  .a  ,F,a 
011223  n-1  n  n 


Proof 

The  path  from  yp  to  yn  implies  that  for  1  <  i  <  n 
the  simulation  state  of  yp  must  be  equal  to  "l”  before  the 
simulation  state  of  yi+1  can  be  set  equal  to  "1",  and  hence 
Fp  must  precede  F^+p  in  any  simulation  sequence.  Q.E.D. 

In  the  example,  the  network  output  graph  contains 
four  paths  of  length  2. 

P(y2,y7)  =  ((y2,y5) (y5,y7)) 
p(yi,,y7)  =  ((y2j,y5)(y5>y7)) 
p(y]|,y7)  = 

FCy^y-^  =  ((yipygMyg^)) 
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By  Property  5 

a 

a 

a 

a 

Every  simulation  sequence  for  the  example  problem 
must  meet  all  four  of  the  above  restrictions.  The  sequences 

ay  =  F2F3F1F2F3 

and 

<xy  =  F2F1F2F3F1 

are  minimum-length  sequences  satisfying  all  four  restrictions. 

Property  5  is  a  necessary  condition  for  ay  of  a  sim¬ 
ulation  sequence.  Property  5  also  provides  a  sufficient  test 
for  an  output  simulation  sequence  if  an  output  simulation  se¬ 
quence  exists.  In  other  words,  if  it  is  known  that  a 
simulation  sequence  exists,  any  sequence  satisfying  Property  5 
is  an  output  simulation  sequence.  There  are  however  networks 
without  a  simulation  sequence,  for  which  an  algorithm  sequence 
satisfying  Property  5  exists. 

Property  6  -  Existence 

An  output  simulation  sequence  a  does  not  exist  for 
a  network  N  if  Ry  is  cyclic. 


=  agF^a-jP^^F3^ 

=  a^F2a^F2agF3cxy 
"  a8F2a9F1a10F3a11 


a12F  0ti3F3a14Fla15 


Proof  -  By  Contradiction 


Assume  Ry  is  cyclic  and  that  a  is  an  output 

simulation  sequence.  Ry  cyclic  implies  that  there  exists 

T  A 

some  y  such  that  yRy.y.  Assume  y  €  F.  Since  s  =  0  and 

V 

s®  =  1,  s  must  be  set  equal  to  1  after  some  F  in  a.  a  is 
thus  of  the  following  form: 


where 


and 


a  = 


Since  yR^y  the  sequence  a-jF  must  contain  F  at  least  twice  and 
is  of  the  form 


a^F  =  a^Fot^F 

where 

ouF 

sy3  -1' 

U-F  a. 

By  Property  s  0  =  1  implies  that  s  =  1.  Q.E.D. 

y  y 

Property  6  states  that  no  output  simulation  sequence 
exists  for  a  network  with  a  cyclic  output  dependency  relation. 
In  the  next  chapter  an  output  simulation  sequence  is  shown  to 
exist  for  any  acyclic  output  dependency  relation. 


IV.  GRAPH  LEVELING 


In  the  previous  chapter  the  problem  of  constructing 
a  simulation  sequence  was  reduced  to  the  problem  of  forming 
an  output  simulation  sequence  for  an  acyclic  output  dependency 
relation.  In  this  chapter  the  structure  of  the  network  out¬ 
put  graph  is  investigated. 

4.1  Leve3  Sets 

Since  simulation  sequences  do  not  exist  for  cyclic 
output  dependency  relations,  the  network  output  graph  is 
assumed  to  be  acyclic.  Any  acyclic  directed  graph  can  be 
graded*  In  this  section  two  types  of  gradings  will  be  de¬ 
fined,  they  are  the  forward  level  sets  and  the  backward 
level  sets. 

Given  any  finite  acyclic  graph  Y  =  <Y,R>  the 
forward  level  sets  L^  are  defined  as  follows: 

L0  *  {y|R(y)  =  $} 
r 

L±  =  •y|R(y)nLi_1  ^  $  and 

i-1 

for  i  >  0  R(y)  C  U  L 

k=0 

The  forward  level  sets  have  a  simple  graphical 
interpretation.  LQ  contains  all  nodes  which  have  no 
predecesscrs.  L^  contains  those  nodes  which  only  have  edges 
q.  A  node  y  is  in  level  L^  (i  >  0)  if  there 


from  nodes  in  L 


to  y  and  every  edge 


exists  an  edge  from  a  node  in  level 
entering  y  comes  from  a  node  in  a  level  less  than  i. 

For  the  example  problem,  the  level  sets  are 

Lo  = 

-  {y^jy^yg^ 

L2  =  {yi»y7> 

L±  =  $  for  i  >  3 

If  Y  is  finite  only  a  finite  number  of  level  sets 
are  nonempty.  Thus  let  SL  equal  the  minimum  i  such  that  a  t. 
Lemma  4.1  -  Only  the  first  i.  forward  level  sets  are  not  empty, 
i.e.  , 

f  *  for  i  <  l 

=  $  for  i  H 

Proof 

By  the  definition  of  £,  ^  $  for  i  <  Si  and 

L.  =  $  for  1  =  1.  Assume  L  =  $  for  some  n  >  i.  L  ...  must 
i  n  —  n+1 

also  be  empty,  since  any  element  of  Ln+^  must  have  a  pre¬ 
decessor  in  L  .  Q.E.D. 

n 

As  a  consequence  of  lemma  ^.1,  £  is  the  number  of 
nonempty  forward  level  sets. 
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Lemma  4.2  -  The  forward-level  sets  are  disjoint,  i.e., 


Li  n  =  $  for  i  /  j 


Proof  -  by  contradiction 

Assume  y  6  O  Lj  for  some  j  >  i.  Since  j  >  0, 
y  €  Lj  implies  that  R(y)Pi  ^j-i  ^  ^  *  Thus  there  exists  at 
least  one  y1  such  that  y1  e  with  (j-1)  ^  i.  But  y  e  L^ 

1  I 

and  y  €  R(y)  implies  that  y  f  for  some  k  <  i.  Q.E.D. 
Lemma  4.3 

The  family  of  sets  for  0  <_  i  <  Jl  partitions  Y. 

Proof 

By  Lemma  4.2,  L^  O  L^  =  $  for  i  ^  j.  Furthermore, 

U  L.  C  Y  by  construction. 
i>0  1 

It  only  remains  to  be  proved  that  Y  C  U  L.  .  Assume 

"  i>0  1 

y  €  Y  and  y  £  VJ  L. .  Clearly  R(y;  f  $>  and  must  contain  some 
i>° 

element  y,  *  (j  L. .  Since  y,  U  L,  ,  R(y,  )  must  contain 

1  i>0  1  1  iVo  1  1 

some  element  y~  t  ■  .  L.  and  this  procedure  may  be  continued. 

*  i^O  1 

Since  Y  is  finite,  there  must  exist  some  y  *  y  such  that 

n  m 


n  n+1  m-1  m 


contradicting  the  acyclic  property  of  R.  Q.E.D. 

The  family  of  sets  L^  thus  partitions  Y.  Each 
element  y  of  Y  is  an  element  of  L^  for  some  i.  The  level  of 
a  nod-' ,  written  Jl(y),  is  defined  as  follows: 
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My)  -  i  if  y  c  ri 


Theorem  4.4 

In  the  system  <Yp,Ry>  with  acyclic  relation  Ry 
aY  =  L0L1L2* 4  *LJl-l 

is  an  output  simulation  sequence  if  each  T.  is  a  block  in  the 
algorithm  partition  of  the  given  network. 

Proof  -  by  induction 

In  order  to  prove  Theorem  4.4  it  is  only  necessary 
to  show  that  each  output  achieves  a  "1"  simulation  state  at 
some  point  in  the  output  algorithm  sequence. 

Consider  any  output  y  such  that  My)  =  0.  This 
output  has  no  predecessors  and  hence  will  achieve  a  "l” 
simulation  state  after  the  algorithm  sequence  a  =  LQ. 

For  all  0  ^  n  £  k  assume  that  Ji(y)  =  n  implies  that 
output  y  will  achieve  a  ”1"  simulation  state  after  the  algorithm 
sequence  a  =  LqL^...L  .  Consider  any  output  y'  such  that 
My')  =  k  +  1.  Any  predecessor  y"  of  y'  is  of  a  level  less 
than  the  level  of  y#.  Thus  My")  <k.  By  inductive 
hypotheses  and  Property  4  the  simulation  statf.  of  /"  is  equal 
to  "1"  after  the  algorithm  sequence  LQL1. .  .L^,.  The  applica¬ 
tion  o  f  Lk+1  will  thus  set  the  simulation  state  of  y '  equal 
to  "1". 


-  50  - 


By  Lemma  4.3  each  output  is  included  in  some 
forward  level  set.  Each  output  will  thus  achieve  a  "1” 
simulation  state.  Q.E.D. 

Theorem  4.4  states  that  the  sequence  of  forward 
level  sets  forms  an  output  simulation  sequence.  This  seauence 
is  in  fact  a  minimum-length  simulation  sequence . 

Property  7  -  Minimum  Length 

The  length  of  any  output  simulation  sequence  ay 
is  greater  than  or  equal  to  l. 

Proof 

To  prove  Property  7,  a  path  P(yQ,yjl_1)  of  length 
8.-1  from  a  node  yQ  f  LQ  to  a  node  y£-1  €  L£_1  will  be 
constructed.  This  path  includes  H  output  nodes.  By 
Property  5,  must  contain  at  least  l  algorithms  and 
hence  is  of  length  greater  than  or  equal  to  l. 

Consider  y£-1  e  L£-i,  by  the  definition  of  the 
level  sets,  there  exists  y£_2  c  L£_2  such  that  y £_2Ryy £_r 
Repeating  the  above  argument,  there  exists  y.  ,  e  L0 
such  that  yj__3^Y^8,-2 *  e4c.  Consider  the  following  sequence: 

y0RYylRY*  *  ,RYyJl-3RYyi.-2RYyjl-l 


This  sequence  is  a  path  containing  l  nodes. 


Q.E.D. 
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As  a  consequence  of  Property  i ,  the  sequence  of 
forward  le/el  sets  forms  a  minimum-length  simulation  sequence. 

It  should  oe  noted  however  that  Theorem  4.4  only  applies  to 
the  special  case  in  which  the  output  algorithm  partition  and 
the  forward  level  partition  are  identical.  This  is  not  the 
case  in  the  example  problem. 

If  an  algorithm  partition  is  desired  for  defining 
the  processes  of  a  network,  the  forward  level  partition  is  an 
excellent  choice.  This  partition  however  is  not  the  only  par¬ 
tition  that  will  generate  a  minimum-length  output  simulation 
sequence.  The  following  method  will  generate  a  partition  which 
is  different  from  the  previous  partition  if  two  or  more  minimum- 
length  output  simulation  sequences  exist. 

Given  a  relation  R  =  {(y,z)}  the  converse  of  R, 

c  c 

written  R  ,  is  the  relation  R  =  {(s,y)>.  If  the  relation 

R  is  represented  by  a  directed  graph  Y,  the  converse  of  R  is 

n 

represented  by  a  directed  graph  Y  obtained  from  Y  by 
reversing  the  direction  of  all  branches. 

Lemma  4 . 5 

0 

R  is  acyclic  iff  R  is  acyclic. 


Proof 

Since  (RC)^  =  R,  it  is  sufficient  to  prove  that 
0 

R  acyclic  implies  that  R  is  acyclic.  Assume  R  Is  acyclic 

c  c 

and  cyclic.  R  cyclic  implies  there  exists  a  loop 
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y1RCy2RC...RCynR°y1 

The  above  implies  that  there  exists 

ylRynR*  *  *Ry2Ryl 

contradicting  that  R  is  acyclic.  Q.E.D. 

Given  a  finite  Y  ^nd  an  acyclic  relation  R  C  Y*Y 

r 

the  backward  level  sets  are  defined  as  follows: 

Lq  =  ( y | RC ( y )  »  *} 

L?  =  ■  y | RC(y)  n  L,  ,  ^  and  RC(y)  C  U 
^  xx  k=0  ^ 

>  . 

c  c 

l  =  minimum  i  such  that  =  $ 

The  reader  will  note  that  the  backward  level  sets 

also  partition  Y.  The  rank  of  the  backward-level  set  parti- 
r 

tion,  l  ,  is  equal  to  Jt,  the  rank  of  the  forward-level  set 
partition. 

The  backward  level  sets  have  the  following  graphical 

r 

interpretation.  A  node  is  in  LQ  if  it  has  no  edges  leaving  it 

r 

A  node  y  is  in  Li  if  there  exists  an  edge  leaving  y  to  a  node 

C 

in  and  every  edge  leaving  y  goes  to  a  node  in  a  level 

less  than  i. 
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For  the  example  problem,  the  backward  level  sets 


Lo  “ 


Li  =  {y3»y5»y6} 


l2  *  ■fy2»y4} 


In  this  example,  the  forward  level  sets  and  the  backward 
level  sets  generate  the  same  partition.  This  is  however 
not  true  in  general. 


Theorem  4.6 


In  the  system  <Yp,Ry>  with  acyclic  relation  Ry 


„  _  Tc  , c  TCTC 

aY  "  hZ-lH-2"  ,L1L0 

is  an  output  simulation  sequence  if  each  Is  a  block  In 

the  algorithm  partition  of  the  given  network. 

Proof 

n 

Consider  an  output  y  €  This  output  has  no 

predecessors.  The  simulation  state  of  y  is  thus  equal  to 

c 

"1"  after  the  algorithm  sequence  a  * 

For  all  1  £  n  £  k  assume  that  y  €  Implies 

that  the  simulation  state  of  y  is  equal  to  ”1"  after  the 

c  c  c 

algorithm  sequence  o  =  L£_1Lj£_2...L£  . 
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P 

Consider  any  output  y  e  LA_^.  Any  predecessor  of 

y  is  in  a  level  greater  than  Z  -  k  and  will  have  a  "1"  simu- 

c  c  c 

lation  state  after  the  algorithm  sequence  a  = 

n 

The  application  of  will  thus  set  the  simulation  state 

of  y  equal  to  Ml".  Q.E.D. 

Corollary  4.7 

The  algorithm  sequence 

M  TC  TC  TCTC 
“y  -  “  ‘"i  0 

is  a  minimum— length  output  simulation  sequence. 

The  application  of  the  backward  level  sets  to  the 
formation  of  output  simulation  sequences  is  also  contingent 
on  the  assumption  that  the  algorithm  partition  and  the  con¬ 
verse  level  partition  coincide.  In  the  next  section  output 
simulation  sequences  will  be  constructed  from  level  partitions 
for  arbitrary  algorithm  partitions. 

4.2  Application  of  Level  Partitions 

In  the  previous  section,  two  partitions,  which 
form  minimum-length  output  simulation  sequences,  were 
developed.  In  this  section  these  partitions  are  applied  to 
generate  output  simulation  sequences  for  arbitrary  algorithm 
partitions.  Before  proceeding  farther  two  substitution  pro¬ 
perties  of  output  simulation  sequences  are  developed. 
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Property  8  -  Substitution 

If  s  aiFa2  is  an  outPut  simulation  sequence  and 
the  output  algorithm  P  C  then  a  =  is  an  output 

simulation  sequence. 

Proof 

By  Property  ^  the  simulation  state  of  an  output 
cannot  change  from  =  1  to  /  1.  The  addition  of  output 
functions  into  a  simulation  sequence  thus  cannot  change  a 
simulation  state  from  =  1  to  ^  1.  Although  this  addition 
may  cause  some  output  simulation  states  to  change  from  ?  1 
to  =  1  at  an  earlier  point  in  the  simulation  sequence,  all 
simulation  states  must  be  *  1  at  the  end  of  the  sequence. 


-  Substitution 


If  oty  =  aiFct2  is  501  outPUt  simulation  sequence 
and  if  output  algorithm  F  =  F2  then  a  *  a^PiF2a2  is  811 

output  simulation  sequence. 


Proof 

Consider  an  output  y  €  FHp^  which  has  "l"  simula¬ 
tion  state  after  the  algorithm  sequence  a^F.  Any  predecessor 
of  this  output  has  ”1"  simulation  state  after  a^.  Output  y 
will  thus  have  "1”  simulation  state  after  a^F^. 

Next  consider  an  output  y  €  F  O  F2  which  has  "1" 
simulation  state  after  the  algorithm  sequence  ct^F.  All 
predecessors  of  y  have  "l"  simulation  state  after  ct-^.  By 
Property  all  predecessors  of  y  have  "1”  simulation  state 
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after  a-^F^ .  Output  y  will  thus  have  "1"  simulation  state 

after  a2FlF2*  Q.E.D. 

In  the  application  of  level  partitions  to  the 
formation  of  an  output  simulation  sequence,  consider  the 
following  three  uses: 

Case  1 

A  level  partition  is  a  refinement  of  the  algorithm 

partition. 

Case  2 

The  algorithm  partition  is  a  refinement  of  a  level 

partition. 

Case  3 

General  case . 

Case  1 

A  partition  L  =  of  a  set  Yp 

is  a  refinement  of  a  partition  Ay  =  {FQ,F1,. . . ,Fn>  of  Yp  if 
for  each  0  £  i  <  £  the  set  is  a  subset  of  Fj  for  some 
0  £  J  n.  Each  set  is  thus  covered  by  some  F^ . 

Theorem  *1.8 

Given  an  output  algorithm  partition  Ay  and  a  forward 
level  partition  L  =  {Lq,L^,. . . if  L  is  a  refinement  of 
Ay  such  that  C  Fi  and  F.^  €  Ay  then  the  algorithm  sequence 
a  =  foP1***^£-1  3-3  a  miniraum“lenSth  output  simulation 


sequence . 
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Proof 

By  Property  8,  a  is  an  output  simulation  sequence. 
This  sequence  is  of  length  Ji  and  by  Property  7  a  minimum- 
length  output  simulation  sequence.  Q.E.D. 

A  special  case  of  Theorem  4.8  is  the  case  of  a 
single  output  algorithm.  Although  this  algorithm  calculates 
new  values  for  every  output,  l  applications  are  required  to 
calculate  correct  values  for  each  output. 

The  reader  will  note  that  a  similar  argument  can  be 
made  for  the  backward-level  partition.  In  fact,  any  partition 
which  forms  an  output  simulation  sequence  can  be  employed. 

Case  2 

If  the  algorithm  partition  Ay  is  a  refinement  of  a 
forward-level  partition  L,  each  output  algorithm  Fi  is  a 
subset  of  some  level  set  L., .  In  this  case  each  level  set 
can  be  expressed  as  the  union  of  cutput  algorithms,  l.e.. 


m(j ) 

L,  =  U  F 


j 


k=l 


j,k  where  PJ ,k  *  pj,i  for  k  *  l' 


Theorem  4.9 


Given  an  output  algorithm  partition  Ay  and  a 
forward  level  partition  L  =  {Lq,L^,.  . .  ,L^^}  if  Ay  is  a 
refinement  of  L  such  that  for  each  0  <  1  <  Z 


L 


i 


m(i ) 
U 
k=l 


9 


where  F±  k  =*  F±  ^  for  k  *  j,  then  the  sequence 


i 

£ 


i 
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°Y  “  F0,lF0,2,**F0,m(0)Fl,l,**Fl,m(l)**'FA-l,m(A-l) 

is  a  minimum-length  output  simulation  sequence. 

Proof 

By  Property  9,  the  sequence  oY  is  an  output  simula¬ 
tion  sequence.  In  this  sequence  each  algorithm  is  included 
exactly  once.  By  Property  2,  the  sequence  oiy  is  a  minimum- 
length  output  simulation  sequence.  Q.E-Tj. 

A  special  case  of  Theorem  ^.9  is  the  case  of  a 
separate  algorithm  for  each  output.  In  this  case  each 
algorithm  is  executed  exactly  once.  The  order  of  execution 
of  algorithms  corresponding  to  outputs  in  the  same  level  set 
is  immaterial.  If  a  multiprocessor  computer  is  available, 
algorithms  corresponding  to  outputs  in  the  same  level  can 
be  executed  concurrently. 

Case  3 

Given  that  Ry  Is  acyclic,  a  simulation  sequence 
can  be  constructed  for  any  algorithm  partition  as  follows: 

1.  Calculate  a  level  partition 

2.  Cover  each  level  set  with  algorithms. 

By  Properties  8  and  9,  the  resulting  sequence  is  an  output 
simulation  sequence . 

Consider  the  example  problem.  For  this  problem 


the  level  sets  are: 


is  an  output  simulation  sequence.  The  above  sequence  contains 
seven  algorithms  and  is  not  of  minimum  length  since  a  sequence 
of  length  5  is  known  to  exist  for  the  example  problem. 
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Although  the  foregoing  method  does  not  yield  minimum- 
length  simulation  sequences,  the  sequences  generated  are 

guaranteed  to  be  simulation  sequences  and  the  effort  required 
is  minimal. 

In  this  chapter  the  graphical  structure  of  the  net¬ 
work  output  graph  has  been  examined.  The  forward  and  backward 
level  partitions  of  the  set  of  outputs  have  been  shown  to 
determine  output  simulation  sequences.  These  partitions 
form  minimum-length  simulation  sequences  for  algorithm 
partitions,  which  are  refinements  of  or  refined  by  a  level 
partition.  The  general  problem  of  forming  a  minimum-length 
simulation  sequence  will  be  discussed  in  the  next  chapter. 
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V.  ALGORITHM  ORDERING 

In  the  previous  chapter  two  partial  orderings  were 
obtained  for  the  set  of  network  outputs .  These  orderings 
specify  output  simulation  sequences  for  any  acyclic  network 
and  minimum-length  output  simulation  sequences  for  a  special 
class  of  algorithm  partitions.  In  this  chapter  an  algorithm 
dependency  relation  is  obtained  from  the  output  dependency 
relation  via  a  condensation  with  respect  to  the  algorithm 
partition. 

In  Section  5.1  an  acyclic  algorithm  dependency 
relation  is  shown  to  specify  a  minimum-length  output  simula¬ 
tion  sequence.  Unlike  the  output  dependency  relation,  the 
algorithm  dependency  need  not  be  acyclic.  A  partitioning  of 
the  ordering  problem  is  discussed  In  Section  5.2. 

5.1  Algorithm  Graph 

The  R^  relation  specifies  tjie  dependency  of  output 
values  on  output  values.  If  y.^  €  F^,  y2  €  F2,  and  y1Ryy2 
the  algorithm  F^  must  precede  the  algorithm  F2  in  any  output 
simulation  sequence.  The  output  dependency  relation  Ry  in 
conjunction  with  the  algorithm  partition  thus  specifies  a 
precedence  relation  among  algorithms.  This  algorithm  relation 
is  formally  defined  as  follows: 


Rp  -  {(F1,F2)  jy1RYy2  for  50me  yl  €  P1  811(1  SOme  y2  €  P2* 
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The  relation  Rp  is  the  condensation  of  Ry  with 
respect  to  the  algorithm  partition.  This  definition  of  con¬ 
densation  differs  from  the  definition  employed  by  Harary  [12] 
in  that  reflexive  elements  are  allowed. 

The  Rp  relation  is  a  precedence  relation  among 
algorithms.  This  relation  is  displayed  by  the  algorithm 
graph  <Ay,Rp>.  A  node  of  the  algorithm  graph  is  an  output 
algorithm,  that  is,  a  nonempty  set  of  outputs.  The  graph  con¬ 
tains  a  branch  from  Fy  to  F2  if  there  exists  at  least  one 
output  y2  c  F2  whose  value  directly  depends  upon  some  output 

yi  €  V 

For  the  example  problem  discussed  in  the  previous 

chapter 

Ay  =  {F1,?2,?3} 

f1  *  {y1,y2,y3> 
f2  *  {yij»y5> 

and 

F3  *  {y6,y7> 

The  network  output  relation  was  shown  in  Figure  5.  The 
algorithm  precedence  relation  Rp  is 

Rp:  (F1^2)  ,  (F1,F3)  ,  (F2,F2)  ,  (F2,F1)  ,  (F2,F3)  ,  (F3,F'L) 


The  algorithm  graph  for  the  example  problem  is  shown 
in  Figure  6a.  The  branches  of  the  network  output  graph  deter¬ 
mining  the  branches  of  the  algorithm  partition  graph  are  shown 
in  Figure  6b.  The  algorithm  partition  graph  is  formed  from 
the  network  output  graph  by  merging  outputs. 

The  reader  will  note  that  Rp,  unlike  Ry,  may  be 
cyclic.  For  the  example  problem  Ry  is  acyclic  and  Rp  is 
cyclic.  If  Rp  is  acyclic,  a  level  partition  of  Ay  can  be 

formed . 

Theorem  5.1 

If  Rp  is  acyclic,  a  forward  level  partition  of  Ay 
determines  a  minimum-length  output  simulation  sequence. 

Proof  (by  induction) 

1.  Consider  an  algorithm  F  €  Lq.  This  algorithm  has 
no  predecessor  algorithms.  Any  output  y  €F  has  no 
predecessor  outputs.  The  application  of  F  will  thus 
set  the  simulation  state  of  all  outputs  in  F  equal 
to  "1."  This  argument  can  be  repeated  for  any 
additional  algorithms  in  LQ.  Thus  if  aQ  is  an 
arbitrary  ordering  of  the  algorithms  in  LQ,  the 
simulation  state  of  every  output  in  an  algorithm 
contained  in  ctg  is  equal  to  "l”  after  Og. 

2.  Let  denote  an  arbitrary  ordering  of  the  algorithms 
in  L*.  Assume  that  all  outputs  in  an  algorithm 
contained  in  the  algorithm  sequence  aoai***an-i 


mitoimaai 
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have  ”1"  simulation  state.  Consider  an  output  y  in 
an  algorithm  F  in  Ln*  Output  y  only  depends  upon 
outputs  included  in  algorithms  in  a  level  less  than 
n.  All  of  these  outputs  have  "1"  simulation  state 
by  inductive  hypothesis.  The  application  of  F  will 
thus  set  the  simulation  state  of  y  equal  to  ’’l”. 

This  argument  can  be  repeated  for  any  additional 
algorithms  in  Ln«  All  outputs  in  an  algorithm  in 
the  algorithm  sequence  aga^...a  thus  will  achieve 
"1"  simulation  state  after  the  application  of  this 
algorithm  sequence. 

In  the  level  partition  each  algorithm  is  In  one 
level  and  hence  occurs  once  in  the  algorithm  sequence.  Since 
each  output  is  in  an  algorithm,  every  output  will  have  "1” 
simulation  state  after  the  algorithm  sequence  Oq0^...o^_^. 

This  sequence  Is  thus  an  output  simulation  sequence  and  by 
Property  2  a  minimum-length  output  simulation  sequence.  Q.E.D. 

The  network  output  graph  displayed  in  Figure  7  has 
an  acyclic  algorithm  dependency  relation.  The  sequences 
F1F2F^F^  and  F^F^F2F^  are  by  Theorem  5*1  minimum-length  out¬ 
put  simulation  sequences. 

The  reader  will  note  that  the  backward  level  sets 
of  Ay  will  also  form  a  minimum-length  simulation  sequence  if 
Rp  is  acyclic. 

If  Rp  is  cyclic,  no  output  simulation  sequence  of 
length  equal  to  the  rank  of  the  algorithm  partitions  will 
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exist  since  at  least  one  algorithm  must  be  applied  more  than 
once . 

5.2  Algorithm  Graph  Partitioning 

Given  any  cyclic  graph  it  is  always  possible  to 
partition  this  graph  into  maximal- connected  subsets.  In  this 
section  the  maximal- connected  subset  partition  is  shown  to 
partition  the  problem  of  forming  a  minimum-length  simulation 
sequence. 

A  set  of  algorithms  K  is  connected  if  for  each 
algorithm  P  in  the  set  there  exist  paths  in  <Ay,Rp>  from  P  to 
every  other  algorithm  of  the  set.  An  empty  set  of  algorithms 
is  not  connected.  A  single  element  set  is  always  connected. 

In  the  algorithm  graph  of  Figure  8  the  following 
sets  are  connected: 


Kx  =  {F^F^F4} 
k2  =  {P1,?2} 

K,  «  {F5} 
k4  =  {F3,F6} 

K5  =  {F1} 

A  maximal  connected  subset  is  a  connected  subset 


which  is  not  contained  in  any  other  connected  subset.  The 
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sets  K^,  ,  and  are  maximal-connected  subsets.  The  set 

K2  is  not  maximal  since  it  is  contained  in  K^. 

Lemma  5.2 

The  family  of  maximal- connected  subsets  is  a 
partition  of  the  set  of  algorithms. 

Proof 

Each  algorithm  is  a  connected  subset  and  hence 
included  in  a  maximal- connected  subset.  The  maximal  con¬ 
nected  subsets  will  be  shown  to  be  disjoint  by  contradiction. 
Assume  and  K2  are  maximal-connected  subsets  with 

F1  cKin  K2  and  F2  €  Ki  -  K2*  Claim  k3  =  UK2  1s  con“ 

nected.  Since  F2  €  K.^  there  exists  a  path  from  F2  to  F^. 

F1  C  K2  iraP1:5-es  fchat  there  exists  a  path  from  F^  to  any 

other  algorithm  in  K2«  Thus  there  exists  a  path  from  F2  to 

any  algorithm  in  K2.  Similarly  there  exists  a  path  from 

any  algorithm  in  K2  to  Fg.  Thus  is  connected.  Q.E.D. 

The  relation  among  maximal  connected  subsets  is 

defined  as  follows: 

iff  i  ^  j  and  F^F^ 

for  some  F^  €  and  some  F^  €  Kj .  is  the  irreflexive 
part  of  the  condensation  of  Rp  with  respect  to  the  maximal 
connected  subset  partition. 

For  the  algorithm  graph  of  Figure  8  the  family  B  of 
maximal  connected  subsets  is 
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B  =  {K1,K3,K2|> 

and  the  RK  relation  among  elements  of  B  is 

V*  (K^),  (K1,K1|),  (K3,K4) 

This  maximal  connected  subset  relation  is  displayed  in 
Figure  9. 

Lemma  5 . 3 

The  relation  Rj^  is  acyclic . 

Proof 

The  reader  will  note  that  Rj^  is  irreflexive  by 
definition.  Assume  that  RR  is  cyclic  and  that  a  cycle  in¬ 
cludes  the  maximal  connected  subsets  K^Kg****,^  for  n  >  1. 
Each  is  connected.  Thus  there  exists  a  path  between  any 
pair  of  algorithms  in  each  connected  subset.  K^Rj^K^,  implies 
that  there  exists  a  path  from  every  algorithm  in  to  every 
algorithm  in  Kg.  The  cycle  thus  implies  that  the  union 
of  the  connected  subsets  is  a  connected  subset.  Q.S.D. 

Since  Rg  is  acyclic,  a  level  ordering  can  be  obtained 
for  the  maximal  connected  subsets. 

Theorem  5.^ 

Given  minimum-length  simulation  sequences  for  each 
maximal  connected  subset  of  A^,  a  forward  level  partition  of 
the  family  of  maximal  connected  subsets  determines  a  minimum- 
length  simulation  sequence. 


•  ' 

(3 


IRE  9 

SUBSET  RELATION 
OF  FIGURE  8 


The  Rk  relation  is  a  precedence  relation  among 
maximal  connected  subsets.  The  proof  of  Theorem  5.4  is 
similar  to  the  proof  of  Theorem  5-1  and  left  to  the  reader. 
Theorem  5*4  partitions  the  problem  of  forming  a  minimum- 
length  simulation  sequence  into  a  collection  of  smaller 
problems . 

Consider  the  network  output  graph  of  Figure  8.  The 

subgraphs  corresponding  to  each  maximal-connected  subset  are 

displayed  in  Figure  10.  The  reader  will  note  that  the 
2  4  12 

sequence  F  F  F  F  is  a  minimum-length  output  simulation  se- 
* 

c 

quence  for  K.^,  that  F^  is  a  minimum-length  sequence  for 
and  that  F^F0  is  a  minimum-length  sequence  for  K^.  As  a 
consequence  of  Theorem  5.4,  F^F^F^F^F^F^F^  is  a  minimum-length 
output  simulation  sequence  for  the  network  displayed  in 
Figure  8. 

In  this  chapter  an  algorithm  dependency  relation 
was  obtained  via  a  condensation  with  respect  to  the  algorithm 
partition.  If  the  algorithm  dependency  relation  is  acyclic, 
a  forward  level  partition  of  the  algorithm  graph  will  deter¬ 
mine  a  minimum-length  simulation  sequence  and  if  the  algorithm 
graph  is  cyclic  the  maximal  connected  subset  partition  will 
split  the  ordering  problem  into  a  collection  of  smaller 
problems. 


NETWORK  OUTPUT  SUBGRAPH  FOR  OUTPUTS 
INCLUDED  IN  ALGORITHMS  IN  K( 


NETWORK  OUTPUT  SUBGRAPH  FOR  OUTPUTS 
INCLUDED  IK  ALGORITHMS  IN  K3 


NETWORK  OUTPUT  SUBGRAPH  FOR  OUTPUTS 
INCLUDED  IN  ALGORITHMS  IN 


FIGURE  10 

NETWORK  OUTPUT  SUBGRAPHS 
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VI .  ALGORITHM  SPLITTING 

In  the  previous  chapter  an  algorithm  ordering 
relation  Rp  has  been  defined.  If  Rp  is  acyclic  a  forward 
level  partition  of  the  algorithm  graph  will  determine  a 
minimum-length  simulation  sequence.  If  Rp  is  cyclic  at  least 
one  algorithm  must  be  included  more  than  once  in  the  output 
simulation  sequence.  In  this  case  some  outputs  of  the 
algorithm  will  achieve  a  ”1”  simulation  state  in  the  first 
application  and  other  outputs  in  the  algorithm  will  achieve 
a  "l"  simulation  state  by  subsequent  applications  of  this 
algorithm.  If  an  algorithm  is  applied  more  than  once  in  a 
simulation  sequence,  this  algorithm  is  equivalent  to  a  sat  of 
algorithms.  One  algorithm  in  this  set  includes  those  outputs 
which  achieve  a  "l"  simulation  state  by  the  first  application 
of  the  algorithm,  another  algorithm  includes  those  outputs 
which  achieve  a  "1"  simulation  state  by  the  second  application 
of  the  algorithm,  etc.  Thus  if  it  is  known  that  an  algorithm 
must  be  included  n  times  In  a  simulation  sequence,  this 
algorithm  can  be  split  Into  n  subsets.  Such  a  split  is 
developed  in  this  chapter. 

In  Section  6.1  output  simulation  sequences  will  be 
characterized  as  refinements  of  algorithm  partitions.  In 
Section  6.2  a  criteria  for  partitioning  an  algorithm  is 
developed.  The  algorithm  partitions  developed  in  Section  6.2 
are  not  unique  and  a  representation  for  this  family  of  level 
partitions  is  presented  in  Section  6.3«  Refinements  of  the 
representation  are  developed  in  Section  6.4. 
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6.1  Output  Partitions 

In  this  section  output  simulation  sequences  will  be 
characterized  as  partitions  of  Yp.  Given  any  minimum-length 
output  simulation  sequence 


let 


“Y  "  P1P2F3- 


and  for 


1  <  i  <  n 


let 


P1P2 


V. 


.  P1F2*“Fi-l  .  , 
and  s  A  pi 

y 


Yi  is  the  set  of  outputs  whose  simulation  state  is 
set  equal  to  ”1"  by  the  application  of  the  ith  algorithm. 

For  the  example  problem  of  Section  2.3 

ay  =  F2F3F1F2F3 

is  an  output  simulation  sequence.  For  this  sequence 

Yi  -  W 


Y2  ’  V 
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y3  =  {yi»y2»y3} 

Yii  =  {y5> 


In  this  section,  the  family  of  Yj^  sets  is  shown  to  be  an 
output  algorithm  partition  refinement  which  forms  an  output 
simulation  sequence. 

Theorem  6.1 

The  family  of  sets  Y.^  for  1  £  i  <_  n  is  a  refinement 
of  the  output  algorithm  partition  Ay. 

Proof 

By  definition,  Y.^  is  the  set  of  outputs  whose 
simulation  state  is  set  equal  to  "1"  by  the  application  of 
P^.  Thus  Yi  Q  for  each  1  <  i  £  n. 

Since  Oy  is  a  simulation  sequence,  the  simulation 
state  of  each  output  y  €  Yp  must  change  from  not  equal  to 
"1"  to  equal  to  "1”  after  the  application  of  some  algorithm. 
Each  output  is  thus  in  an  element  of  Y^  for  some  1  <  i  £  n. 

Q .  £ .  D . 


Theorem  6.  2 

The  sequence 


a  =  Y.Y9. . . Y 
12  n 

is  an  output  simulation  sequence. 


J 


V 
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Proof 

Since  the  family  Y^  partitions  Yp,  each  output 
y  €•  Yp  is  an  element  of  Y^  for  some  j  and  will  achieve  a 
"1"  simulation  state  after  the  application  of  Yj .  By 
Property  4,  this  output  will  have  "l”  simulation  state  after 
the  sequence  Y^Yg^.Y  .  Q.E.D. 

Given  any  minimum-length  output  simulation  sequence 
of  length  n,  a  partition  of  Yp  of  rank  n  can  be  extracted  from 
the  simulation  sequence  such  that  the  partition  will  form  a 
simulation  sequence.  In  this  simulation  sequence  no  algorithm 
is  executed  more  than  once.  Thus  for  every  network,  for  which 
a  simulation  sequence  exists  (Ry  is  acyclic),  there  exists  a 
partition  which  is  a  refinement  of  Ay  and  which  can  be  ordered 
to  form  a  minimum-length  output  simulation  sequence.  The 
problem  of  finding  a  minimum-length  output  simulation  sequence 
is  thus  equivalent  to  the  problem  of  finding  a  minimum-rank 
partition  which  is  a  refinement  of  Ay  and  for  which  Rp  is 
acyclic . 

6.2  Compatible  Output  Sets 

In  this  section  a  criteria  for  partitioning  algorithms 
is  developed.  Suppose  that  an  algorithm  contains  two  Ry  re¬ 
lated  outputs.  By  the  path  ordering  property,  this  algorithm 
must  be  included  at  least  twice  in  any  output  simulation 


sequence . 
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In  the  example  problem  of  Section  2.3  outputs  y^ 

? 

and  y<-  are  R^-related  and  are  elements  of  algorithm  F  .  If  the 
2 

algorithm  P  Is  split  into  two  algorithms  and  where 

*  p2  -  =  ^y5> 

P2  =  F2  -  {y5)  =  {y^> 

and  a  minimum-length  simulation  sequence  is  obtained  for  the 
new  algorithm  partition,  this  sequence  will  be  a  minimum- 
length  simulation  sequence. 

In  general  there  can  be  many  pairs  of  -related 
outputs .  A  compatible  output  set  is  a  set  of  outputs  con- 
taining  no  pair  of  Ry-related  outputs.  A  singleton  output 
set  is  a  compatible  output  set  for  any  acyclic  output  de¬ 
pendency  relation.  The  reader  will  note  that  any  partition 
representing  a  minimum-length  output  simulation  sequence 
will  contain  only  compatible  sets.  Thus  if  an  algorithm 
is  not  a  compatible  set,  only  subsets  of  this  algorithm 
need  be  considered  in  the  formation  of  an  output  simulation 
sequence. 

For  each  output  algorithm  F  consider  the  output 
+ 

algorithm  graph  <F,Ry>.  This  graph  represents  the  total 
dependency  relation  between  outputs  in  the  algorithm  F. 

The  output  algorithm  graphs  for  the  three  algorithms 
of  the  example  problem  are  shown  in  Figures  11a,  lib,  and  11c. 

p  4 

Only  F  contains  Ry -related  outputs. 
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*  yl 

<f\ry>  -  y2 

“  y3 


(a)  Output  Algorithm  Graph  for  F1 


yH 


(b)  Output  Algorithm  Graph  for  F” 


<F3,  Rj> 


y6  y7 


(c)  Output  Algorithm  Graph  for  F3 


FIGURE  11 

Output  Algorithm  Graphs  for 
Algorithms  of  Figure  6 
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Theorem  6 . 3 

+ 

The  forward-level  partition  of  <F,Ry>  is  a  minimum 
rank  partition  of  F  into  compatible  sets. 

Proof 

4* 

If  two  outputs  are  Ry  related  they  must  be  in 
separate  levels.  Each  forward-level  set  is  thus  a  compatible 
set.  To  prove  minimum  rank,  suppose  that  the  output  algorithm 
graph  contains  a  maximal -length  path  containing  %  nodes. 
Clearly  any  compatible  set  can  contain  at  most  one  of  these 
nodes.  The  rank  of  the  partition  is  thus  greater  than  or 
equal  to  l.  The  rank  of  the  forward-level  partition  is 
equal  to  I.  Q.E.D. 

In  a  similar  manner  it  can  be  proved  that  the 

+ 

backward  level  set  partition  of  <F,Ry>  is  a  minimum-rank 
partition  of  F  Into  compatible  sets.  In  fact,  the  forward- 
and  backward  level  set  partitions  are  only  two  of  a  family 
of  minimum-rank  level  partitions.  The  choice  of  minimum- 
rank  partition  is  relevant  to  the  construction  of  a  minimum- 
length  simulation  sequence. 

6.3  Maximal  Algorithm  Subset 

In  this  section  a  representation  for  the  family  of 
level  partitions  will  be  developed.  The  number  of  sets  in  a 
level  partition  has  previously  been  shown  to  be  equal  to  the 
number  of  nodes  in  the  longest  path.  By  the  path  ordering 
property , outputs  on  the  same  path  must  be  in  different  level 
sets.  Each  output  of  a  maximal -length  path  is  in  a  separate 
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level  set  and  each  level  set  contains  one  output  of  each 
maximal-length  path. 

For  the  purpose  of  discussion  it  will  be  assumed 
that  the  level  sets  are  numbered  in  the  order  in  which 
they  will  form  a  simulation  sequence .  The  first  forward  level 
set  is  thus  Lq,  and  the  first  backward-level  set  is 

Suppose  that  an  output  is  an  isolated  node,  that  is, 
+ 

it  is  not  Ry  related  to  any  other  output  in  the  algorithm. 
Clearly  this  output  can  be  included  in  any  level  set.  In  a 
forward  level  partition  isolated  outputs  are  included  in  the 
first  level  set.  In  the  backward-level  partition  isolated 
outputs  are  included  in  the  last  level  set .  In  fact ,  the 
forward  level  partition  has  the  property  that  each  output  is 
included  in  the  earliest  level  set  and  the  backward-level 
partition  has  the  property  that  each  output  is  included  in 
the  latest  level  set. 

Consider  the  output  algorithm  graph  of  Figure  12, 
the  forward  and  converse  level  partitions  of  this  graph  are: 


LQ  =  (1,3, 8, 9) 
Lx  =  {2,4} 

L2  =  {6} 

L3  =  {7) 

=  {5> 

Lc  =  {10} 


lP5  =  {3} 

4  =  W 

Lc3  =  {6,8} 

lP2  =  {7} 

=  {1,5> 

Ln  =  {2,9,10} 
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This  output  algorithm  graph  has  a  path  containing 
six  nodes,  i.e.,  the  path 


3+4+6-7-5-10 


Hence  the  algorithm  must  be  included  at  least  six  times  in  any 
output  simulation  sequence.  After  the  first  application  of 
this  algorithm  only  nodes  in  LQ  can  achieve  a  "1"  simulation 
state.  In  the  initial  application  of  this  algorithm  only 
nodes  in  Lq  need  be  considered.  If  this  algorithm  is  to  be 
included  only  six  times  in  the  output  simulation  sequence, 
the  nodes  of  Lq  which  are  on  a  maximal-length  path  must 
achieve  "1"  simulation  state  after  LQ.  These  nodes  are  also 

c 

in  In  the  output  algorithm  graph  of  Figure  12  only 

node  3  must  achieve  a  "1"  simulation  state  after  the  first 
application  of  the  algorithm.  Nodes  1,  8,  and  9  can  achieve 
a  ”1"  simulation  state  after  a  later  application.  These  nodes 
can  achieve  a  "1"  simulation  state  after  the  second  application. 
Additional  nodes  which  can  achieve  a  "l"  simulation  state 
after  the  second  application  are  the  nodes  in  L^.  After  the 
second  application  of  the  algorithm  in  any  sequence  containing 
only  a  minimum  number  of  these  algorithms,  the  set  of  nodes 
which  can  attain  a  ”1"  simulation  state  is  (Lq  - 
Only  nodes  in  the  above  set  need  be  considered  in  the  second 
application. 
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In  general  the  set  of  outputs  which  can  attain  a 
"1"  simulation  state  after  the  nth  application  of  an 
algorithm  in  any  sequence  containing  only  a  minimum  number 
of  these  algorithms  is  called  a  maximal  subset  of  an 
algorithm.  A  maximal  subset  is  not  necessarily  a  compatible 
set.  For  the  output  algorithm  graph  of  Figure  12 

t 

(L0  “  LLl)  UL1  = 

This  set  is  not  a  compatible  set  since  nodes  1  and  2  are 
Ry  related.  Outputs  1  and  2  are  on  a  path  containing  two 
nodes.  Since  at  least  six  applications  of  the  algorithms 
are  required,  any  one  of  the  first  five  applications  may  set 
the  simulation  state  of  output  1  equal  to  "1". 

Given  the  forward  and  backward-level  partitions 
the  family  Cj  of  maximal  subsets  is  formed  as  follows: 


C!  *  (C0  -  L*-l]  U  h 

ci  -  h-i  -  u  l. 
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For  the  acyclic  output  algorithm  graph  of  Figure  12,  the 
maximal  subsets  are: 

C0  =  U, 3,8,9) 

Cx  =  {1,2, 4, 8, 9) 

C2  =  {1,2, 6, 8, 9) 

C3  =  {1,2, 7, 9) 

C4  =  {1,2, 5, 9) 

C5  =  {2,9,10} 

The  family  of  maximal  algorithm  subsets  is  a 
cover  of  the  family  of  minimum-rank  level  partitions,  that  is, 
each  forward  level  set  is  a  subset  of  the  corresponding 

Q 

maximal  subset  and  each  backward  level  set  is  a  subset 
of  Cfl,-i-l‘ 

6.4  Refinement  of  Maximal  Algorithm  Subsets 

In  Section  6.3,  a  family  of  maximal  subsets  was 
defined  for  each  algorithm.  These  subsets  are  a  cover  of 
the  family  of  minimum-rank  level  partitions.  In  the  maximal 
subsets  only  nodes  on  a  maximal -length  path  in  the  output 
algorithm  graph  are  included  in  a  single  unique  subset.  In 
this  section  refinements  of  maximal  algorithm  subsets  are 


developed. 
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An  output  of  an  algorithm  is  ordered  if  it  is 
included  in  exactly  one  maximal  subset.  The  output  algorithms 
graph  of  Figure  12  contains  one  maximal-length  path  with  the 
following  six  nodes: 

3+-4+6->-7  +  5+10 

The  above  nodes  are  ordered.  The  remaining  nodes  are  included 
in  more  than  one  maximal  subset  and  hence  not  ordered. 

In  the  formation  of  maximal  algorithm  subsets  only 
the  dependency  of  outputs  on  outputs  in  the  same  algorithm 
was  considered.  In  this  section  the  structure  of  the  network 
output  graph  is  examined. 

The  predecessor  set  of  an  output  y, written  P(y), 
is  the  set  of  outputs  which  have  a  path  to  output  y. 

P(y)  -  jy' |y'Ryy} 

The  successor  set  of  an  output  y,  written  S(y), is  the  set  of 
outputs  which  are  accessible  from  output  y. 

s(y)  =  {y'iyRyy'} 

Consider  the  network  output  graph  of  Figure  13*  The  output 
algorithm  graph  of  algorithm  F  was  shown  in  Figure  12 .  The 
maximal  algorithm  subsets  were  developed  in  the  previous 
section.  Output  8  is  included  in  the  first  three  subsets. 


NETWORK  OUTPUT  GRAPH 
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Since  it  has  no  predecessors,  the  first  application  of 
algorithm  will  set  its  simulation  state  equal  to  "l". 
Output  8  thus  need  be  included  only  in  Cq. 

Rule  1 

If  the  predecessor  set  of  an  output  is  empty.  In¬ 
clude  this  output  only  in  CQ. 

Next  consider  output  9.  This  output  has  no 
successors.  Thus  no  output  depends  upon  this  output.  In 
this  case  it  is  desirable  to  include  output  9  in  the  last 
maximal  subset. 

Rule  2 

If  the  successor  set  of  an  output  is  empty,  in¬ 
clude  this  output  only  in  the  last  subset. 

Output  1  has  the  same  predecessor  set  as  output  3. 
Since  outputs  1  and  3  are  in  the  same  algorithm,  any  algorithm 
sequence  which  sets  the  simulation  state  of  one  output  equal 
to  ”1"  will  also  set  the  simulation  state  of  the  other  output 
equal  to  "l".  These  outputs  can  thus  be  included  in  the 
same  subset.  Output  3  is  ordered  and  included  only  in  Cq. 
Output  1  can  thus  be  ordered. 

Rule  3 

If  output  y  is  included  only  in  subset  n  and 
P(y)  =  P(y#)  include  output  y#  only  in  subset  n. 

Output  2  and  output  10  have  the  same  successor  set . 
All  outputs  which  are  dependent  upon  output  2  are  also 
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dependent  upon  output  10  and  vice  versa.  Outputs  2  and  10  can 
thus  be  included  in  the  same  subset. 

Buie  4 

If  output  y  is  included  only  in  subset  n  and 
S(y)  =  S(y')  include  output  y'  only  in  subset  n. 

The  ordering  rules  presented  here  are  frequently 
contradictory.  For  example,  rule  2  can  be  applied  to  include 
output  9  into  subset  or  rule  3  can  be  applied  to  include 
output  9  into  subset  Cq.  Either  application  can  lead  to  the 
formation  of  a  minimum-length  output  simulation  sequence. 

For  the  network  output  graph  of  Figure  13,  the 
refined  maximal  algorithm  subsets  are: 

CQ  =  {1,3,8} 

c3  =  {4} 

C2  =  {6} 

C3  =  (7> 

=  (5) 

C5  =  (2,9,10) 

The  network  output  graph  with  the  above  subsets  substituted 
for  F^  is  shown  in  Figure  l1*.  The  algorithm  dependency 
relation  for  this  graph  is  acyclic.  The  sequence 
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p^p3p1p3p3p5p3p3F3F2 

is  a  minimum-length  output  simulation  sequence. 

In  this  example  problem  the  algorithm  splitting 
technique  has  resulted  in  an  acyclic  algorithm  graph. 

This  is  not  always  the  case.  In  the  next  chapter  a  general 
method  for  forming  simulation  sequences  is  presented. 


VII.  FORMATION  OF  SIMULATION  SEQUENCES 

In  previous  chapters  various  techniques  for  con¬ 
structing  simulation  sequences  were  presented.  In  Chapter  III 
the  problem  of  constructing  a  simulation  sequence  was  reduced 
to  the  problem  of  forming  an  output  simulation  sequence.  In 
Chapter  V  an  algorithm  ordering  relation  was  discussed  and 
the  simulation-sequence  formation  problem  was  reduced  to  a 
collection  of  smaller  problems.  A  criterion  for  splitting 
algorithms  was  developed  in  Chapter  VI.  In  this  chapter  these 
techniques  are  combined  into  a  general  method  for  construct¬ 
ing  a  minimum-length- simulation  sequence. 

In  Section  7.1  a  selective  search  technique  is 
presented.  The  general  method  is  presented  in  Section  7.2 
and  an  illustrative  example  is  presented  in  Section  7.3. 

7.1  Modified  Eligible-Set  Method 

In  this  section  a  selective  search  technique  for 
the  construction  of  a  simulation  sequence  will  be  developed. 

Although  this  technique  is  applicable  to  all  acyclic  networks, 
it  is  primarily  intended  for  a  problem  to  which  none  of  the 
techniques  discussed  in  Chapters  V  and  VI  apply,  that  is,  problems 
represented  by  a  connected  algorithm  graph  of  maximal  algo¬ 
rithm  subsets . 

In  Section  2.3  a  simulation  sequence  was  constructed 
via  a  sequential  selection  of  algorithms.  The  choice  of 
algorithms  was  restricted  by  a  one-s^ep  look-ahead  rule. 
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The  eligible  output  set  Ey  Includes  those  outputs  which  do 
not  have  "1"  simulation  state  but  which  can  achieve  ”1" 
simulation  state  via  a  single  application  of  an  algorithm. 
Formally 


ey  =  *yisy  t  i  and  y,Ryy  =>  sy>  s  i) 

The  eligible  algorithm  set  EA  is  the  set  of  algorithms  which 
have  at  least  one  output  in  the  eligible  output  set.  Formally 

EA  =  (F |  y  c  Ey  for  some  y  c  F} 

Both  eligible  sets  are  a  function  of  the  algorithm 
sequence  previously  applied.  In  the  search  for  simulation 
sequences  only  algorithms  in  EA  need  be  considered.  A 
minimum-length  simulation  sequence  can  be  constructed  by 
considering  all  eligible  algorithms  at  each  step.  This 
method  will  require  the  construction  of  a  large  number  of 
algorithm  sequences.  The  computation  time  can  be  reduced 
drastically  by  recognizing  optimal  choices.  Suppose  that 
the  entire  set  of  outputs  of  an  algorithm  is  included  in 
the  eligible  output  set .  Selecting  this  algorithm  as  the 
next  algorithm  in  the  sequence  is  an  optimal  choice  since  the 
algorithm  will  be  included  only  once  in  the  simulation 
sequence.  Each  output  in  the  algorithm  will  achieve  a  "l" 
simulation  state  and  outputs  which  are  dependent  upon  outputs 
in  this  algorithm  may  become  eligible. 


The  done  output  set  Dy  is  the  set  of  outputs  which 
have  "l”  simulation  state.  The  done  algorithm  set  is  the 
set  of  algorithms  which  have  all  their  outputs  in  the  done 
output  set.  A  selection  of  an  algorithm  is  an  optimal  choice 
if  the  select ion^ adds  this  algorithm  to  the  done  algorithm 
set.  The  set  of  maximum  algorithms  includes  those  al¬ 
gorithms  in  Ea  which,  if  selected, will  become  elements  of  DA. 

The  order  of  choice  of  an  algorithm  from  a  set  of  maximum 

algorithms  is  immaterial. 

If  the  eligible  algorithm  set  includes  only  a 
single  algorithm,  the  selection  of  this  algorithm  is  by 
default  an  optimal  choice.  In  the  case  of  two  or  more 
algorithms,  each  choice  must  be  considered.  This  can  be  ac¬ 
complished  by  a  tree  structure  search.  The  flow  chart  dis¬ 
played  In  Figure  15  employs  such  a  technique. 

The  parameter  L  indicates  the  level  of  the  tree 
structure.  Initially  L  is  equal  to  1.  L  is  increased  by 
1  whenever  a  nonoptimal  choice  is  encountered.  L  is  de¬ 
creased  by  1  after  all  choices  in  the  largest  level  have 
been  exausted. 

Before  incrementing  L,  a  test  Is  performed  to 
check  whether  the  current  algorithm  sequence  can  form  a  simu¬ 
lation  sequence  of  length  less  the  shortest  simulation  sequence 
found.  The  potential  length  of  a  simulation  sequence  is  equal 
to  the  number  of  algorithms  in  the  sequence  plus  the  number 


SELECT  ALGORITHM 
IN  Ea  (L-l) 

INITIALIZE  LEVEL  L 


FIGURE  15 

FLOW  CHART  FOR  THE  SELECTIVE  SEARCH  METHOD 


-  96  - 


of  algorithms  which  are  not  in  D^.  If  the  potential  length 
is  greater  than  or  equal  to  MINL,  this  sequence  cannot  form 
a  minimum-length  simulation  sequence.  In  this  case  the  routine 
abandons  this  sequence  and  selects  the  next  choice  in  the 
highest  level.  If  no  choices  remain  and  L  is  equal  to  1  the 
routine  stops.  If  L  is  greater  than  1,  L  is  decreased  by 


Whenever  an  algorithm  is  selected, this  algorithm 
is  removed  from  the  eligible  algorithm  set  and  added  to  the 
simulation  sequence.  The  outputs  in  the  algorithm  which 
achieve  nlM  simulation  state  are  added  to  the  done  output 
set.  A  new  eligible  set  is  formed  and  the  procedure  is 
continued. 

If  the  eligible  output  set  or  the  eligible  algorithm 
set  is  empty,  all  outputs  which  car.  achieve  "1"  simulation 
state  have  achieved  ”1”  simulation  state  and  the  length  of 
the  sequence  formed  is  compared  with  the  shortest  sequence 
found.  The  reader  will  note  that  this  sequence  is  a  simula¬ 
tion  for  any  acyclic  network.  In  the  case  of  a  cyclic  net¬ 
work  the  sequence  formed  will  result  with  a  "l”  simulation 
state  for  all  outputs  which  are  not  dependent  upon  outputs 
in  a  cycle.  The  modified  eligible  set  method  thus  forms  a 
simulation  sequence  for  acyclic  networks  and  in  the  case  of 
cyclic  networks  this  method  will  form  a  sequence  which  will 
result  In  "1”  simulation  state  for  all  outputs  which  can 
achieve  "1”  simulation  state. 
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Consider  the  example  problem  discussed  in  Sec¬ 
tion  2.3.  The  network  output  graph  for  this  problem  was 
shown  in  Figure  5,  the  algorithm  graph  was  given  in  Figure  6, 
and  the  output  algorithm  graphs  were  given  in  Figure  11. 
Algorithm  F  is  split  into  two  algorithms  by  the  algorithm 
splitting  technique.  The  network  output  graph  with  maximal 
subsets  substituted  for  F2  is  shown  in  Figure  16.  In  order 
to  simplify  the  notation  each  node  has  been  labeled  with 
the  output  subscript. 

Initially  L  =  l,  MINL  =  7  (number  of  outputs), 
SEQ(l)  =  \ ,  Dy(l)  =  $  and  DA(1)  =  4>.  The  eligible  sets  are 

Ey  =  ,  and  =  {f2>. 

Since  Em  is  not  empty,  the  algorithm  in  E^  is  selected 

» 

SEQ(l)  =  F2,  Dy (l)  =  (20  and  DA(1)  =  {F2>. 

The  eligible  sets  become: 

Ey  =  {2*3,6},  EM  =  $  and  EA(l)  *  {F1,F^}  . 

A  choice  exists.  The  potential  length,  4 ,  is  less  than  7. 

Thus  consider  each  choice,  set  L  =  2,  select  F1,  and 
initialize  the  new  level. 

SEQ(2)  =  F1,  Dy(2)  =  {2,3,4}  and  Dft(2)  =  (F2). 


FIGURE  16 

NETWORK  OUTPUT  GRAPH  OF  FIGURE  6b 
WITH  MAXIMAL  SUBSETS  SUBSTITUTED  FOR  F2 


-  99  - 


The  eligible  sets  are  now 

Ey  =  {5,6}  and  Em  =  {F.^. 

E^  is  non-empty,  thus  select  F^. 

SEQ(2 )  =  F1F1,  Dy(2)  =  {2, 3, 4,5)  and  DA(2)  =  {F^F^. 

The  new  eligible  sets  are 

Ey  =  {6,7}  and  EM  =  {F3} 

-j 

F  is  a  maximum  algorithm. 

S£Q(2)  =  F1F1F3,  Dy(2)  =  (2, 3, 4, 5, 6,7),  DA(2)  =  {F1,F2,F3}, 

Ey  =  {1},  Em  =  {F1}. 

F1  is  a  maximum  algorithm  since  only  output  1  remains  to  be 
done. 

SEQ(2)  =  F1P1F3F1i  Dy(2)  =  {1 ,2,3,4 ,5,6,7),  DA(2)  =  {F1,F1,F2,F3} 
and  Ey  3  $  . 

1  1 

The  sequence  FgF  F^F-'F  is  a  simulation  sequence. 

MINI  =  5 

The  routine  now  selects  F  ,  initializes  level  2  and  proceeds 


as  follows : 


SEQ(2)  *  F3  Dy(2)  *  {4,6}  DA(2)  =  {F2) 

Ey  =  {1,2,3}  Em  =  {F1} 

SEQ(2 )  =  p¥  Dy(2)  =  {1,2, 3, 4, 6}  DA(2)  *  {F1,F2} 

Ey  =  {5}  Em  =  {Fx} 

SEQ(2)  =  F3F1F1  Dy(2)  =  {1,2, 3, 4, 5,6}  DA(2)  =  {F1,F1,F2} 

Ey  =  {7}  Em  =  {F3} 

SEQ ( 2 )  *  F3F1F1F3  Dy(2)  =  {1,2, 3, 4, 5, 6, 7} 

Da(2)  =  {F1,F1,F2,F3} 

Ey  =  $ 

•3  1  3 

The  new  simulation  sequence  is  FgF  F  F-^F  .  In  this  case  no 
choices  remain  and  the  routine  terminates. 

The  eligible  set  algorithm  can  form  a  minimum-length 
simulation  sequence  for  any  network  with  acyclic  Ry.  The 
algorithm  forms  the  simulation  sequence  by  a  sequential  selec¬ 
tion  of  algorithms.  A  one-step  look  ahead  technique  is 
employed.  Optimal  choices  are  recognized  and  a  potential- 
length  test  abandons  simulation  sequences  whose  length  will 
be  greater  than  the  minimum-length. 
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The  amount  of  computation  that  is  required  is 
dependent  upon  the  problem.  If  the  algorithm  graph  is 
acyclic,  the  maximum  algorithm  set  is  never  empty  and  the 
number  of  eligible  sets  formed  is  equal  to  the  number  of 
nodes  in  the  longest  path  of  the  algorithm  graph.  For  the 

network  of  Figure  7,  the  routine  will  select  F  in  the 

2  ii  q 

first  step,  F  and  F  in  the  second  step,  and  FJ  in  the 

third  and  final  step.  The  routine  is  thus  very  efficient 

for  acyclic  algorithm  graphs. 

In  the  case  of  a  cyclic  algorithm  graph  the  amount 
of  computation  is  dependent  upon  the  choice  of  algorithms. 

The  routine  will  examine  all  potential  minimum-length 
simulation  sequences.  In  the  routine  a  potential  minimum- 
length  simulation  sequence  is  any  simulation  sequence  which 
is  shorter  than  the  shortest  sequence  found.  Although  the 
routine  will  halt  with  a  minimum-length  sequence,  the  amount 
of  computation  Is  dependent  upon  the  length  of  sequences 
found.  If  the  initial  selection  results  in  a  short  sequence, 
a  large  number  of  sequences  will  be  abandoned  by  the  potential 
length  test . 

The  amount  of  computation  can  be  reduced  by  Incor¬ 
porating  simple  selection  rules  which  by  ordering  eligible 
algorithms  will  cause  the  search  technique  to  examine  the 
'  more  likely  sequences  first.  Prior  to  any  selection  of  a 
nonmaximum  algorithm  each  eligible  algorithm  is  examined 
and  the  best  algorithm  is  selected.  The  following  selection 
rules  have  been  investigated  by  the  author: 


102 


1.  Select  the  algorithm  which  will  update  the  most 
outputs . 

2.  Update  the  output  in  the  lowest  forward  level. 

3.  Update  the  output  in  the  highest  backward  level. 

None  of  these  rules  apply  to  all  problems.  The 
author  has  constructed  problems  to  which  none  of  the  above 
rules  apply.  Rule  1  is  based  on  the  presumption  that  if  the 
maximum  number  of  outputs  is  updated  at  each  step ,  the 
number  or  steps  will  be  minimized.  The  presumption  is  false 
and  Rule  1  is  not  a  good  selection  criteria  since  output 
fan-outs  are  not  considered.  The  reader  will  note  that  if 
an  output  does  not  fan  out,  then  whether  this  output  is  up¬ 
dated  at  the  present  step  or  at  a  future  step  is  immaterial. 

Rules  2  and  3  ere  based  on  the  path  criterion,  that 
is,  outputs  that  are  initial  nodes  of  paths  in  the  network 
graph  must  be  updated  before  their  successors  can  be  up¬ 
dated.  The  level  partitions  of  the  network  output  graph 
specify  such  an  ordering.  Rule  3  is  usually  a  better  selec¬ 
tion  rule  than  Rule  2  since  in  Rule  3  outputs  that  are 
initial  nodes  of  maximal-length  paths  are  given  preference. 

The  rules  mentioned  above  are  simple  examples  of 
selection  rules.  More  powerful  selection  rules  can  be 
developed.  The  implementation  of  selection  rules  will  reduce 
*  the  amount  of  computation  only  if  the  computation  required 
by  the  selection  rule  is  less  than  the  computation  required 
by  the  eliminated  sequences.  It  is  the  opinion  of  the 
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author  that  except  for  special  uses  the  implementation  of  selectior 
rules  will  not  significantly  reduce  the  total  amount  of  computation. 

7.2  General  Method 

Although  the  modified  eligible- set  method 
presented  in  the  previous  section  can  form  mlnimum-lgteth 
simulation  sequences  for  any  acyclic  network,  it  doe^Bfct 
employ  the  maximal  connected  subset-  partitioning  teefiapue 
developed  in  Chapter  V.  In  this  section  a  general  meJiPDd 
is  presented.  This  method  is  simply  an  ordering  of  the 
techniques  previously  discussed. 

1.  Calculate  maximum  algorithm  subsets. 

2.  Form  an  algorithm  graph  using  the  maximum  algorithm 

suosets . 

3.  Partition  the  algorithm  graph  into  maximal-eonnected 
subgraphs. 

4.  Level  order  the  maximal- connected  subsets. 

b.  Obtain  a  minimum-length  simulation  sequence  for 
each  maximal- connected  subset. 

In  the  general  method  the  algorithm  splitting 
technique  discussed  in  Chapter  VI  is  applied  first.  This 
technique  results  in  a  new  set  of  algorithms.  The  partition¬ 
ing  of  the  ordering  problem  discussed  in  Chapter  V  is 
applied  next  and  the  modified  eligible  set  method  is  applied 
last  to  each  maximal  connected  subset. 

A  computer  program  has  been  written  in  the 
FORTRAN  V  language.  This  program  forms  simulation  sequences 
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using  the  general  method.  The  program  contains  approximately 
1400  cards  and  is  operational  on  the  Uni vac  1108  computer.* 

The  operation  of  the  program  in  the  solution  of  an  example 
problem  is  described  in  the  next  section. 

7.3  Illustrative  Problem 

In  this  section  a  simulation  sequence  will  be  con- 
sfl|w;ed  for  an  illustrative  problem.  The  system  to  be 
simulated  is  composed  of  12  subsystems.  Each  subsystem 
contains  memory  and  at  least  two  output  variables.  The  pro¬ 
cess  variable  dependency  relation  is  assumed  to  be  known  for 
each  process.  The  network  interconnections  are  also  known. 

The  behavior  of  each  process  is  described  by  two  algorithms, 
one  algorithm  calculates  new  values  for  the  state  variables 
and  a  second  algorithm  calculates  new  values  for  the  output 
variables. 

The  output  dependency  relation  for  this  network 
is  shown  in  Figure  17.  The  outputs  and  processes  have  been 
numbered  numerically.  The  algorithm  definition  is  stored  by 
the  computer  program  in  the  matrix  ALG.  This  matrix  is  shown 
in  Figure  18.  The  output  dependency  relation  is  stored  in  the 
Boolean  matrix  OUT.  The  entry  in  row  i  and  column  j  of 
this  matrix  is  equal  to  1  if  output  j  directly  depends  upon 
output  i.  The  OUT  matrix  is  shown  in  Figure  19-  The  row  2, 
column  5  entry  of  this  matrix  is  1  since  output  5  depends 
upon  output  2. 
l - 

A  listing  of  c-n.'  program  is  available  from  the  author.  A 

user  bul^tin  for  this  program  is  presented  in  the  appendix. 
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ROW  NO.  ROW  ENTRIES  MATRIX  ALG 


I 

1 

2 

3 

2 

5 

6 

3 

7 

8 

4 

9 

10 

5 

11 

12 

13 

6 

14 

15 

16 

7 

17 

18 

19 

8 

20 

21 

22 

9 

30 

31 

10 

28 

29 

11 

26 

27 

12 

23 

24 

25 

FIGURE  18 

Algorithm  Definition 
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MATRIX  OUT 


1  2  3 

ROW  1 2345678901 2345678901 2345678901 

1  0000000000000000000000000000000 

2  00001 00000000000000000000000000 

3  000001 0000000000000000000000000 

4  0000000000000000000000000000000 

5  0000001 000000000000000000000000 

6  00000001 00000000000000000000000 

7  000000001 0000000000000000000000 

8  0000000001 000000000000000000000 

9  0000000000000000000001 000000000 

1 0  00000000000000000001 00000000000 

1  1 000000000000000000000000000000 

2  0000000000000000000000000000000 

3  0000000000000000000000000000001 

4  0000000000001 000000000000000000 

5  00000001 00000000000000000000000 

6  0001 000000000000000000000001 000 

7  0000000000000001 000000000000000 

8  0000000001 000000000000000000000 

9  000000000000000000000001 01 00000 

20  0000000000000000001 000000000000 

1  00000000000000000000000001 00000 

2  00000000000000000000001 00000000 

3  0000000000000000000000000000000 

4  OOOOOOOOOOOCOOOOOOOOl 0000000000 

5  000000000000000000000000001 0000 

6  00000000000000000000000000001 00 

7  0000000000000000000000000000000 

8  0000000000000100000000000000000 

9  00000000000000001 00000000000001 

30  000000000001 0000000000000000000 

1  0000000000000000000000000000000 


FIGURE  19 

Output  Dependency  Relation  Matrix 
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Each  entry  in  this  matrix  represents  one  bit  of  a 
computer  memory  word.  Since  the  Univac  1108  computer  has 
36  bits  per  word,  each  row  of  the  matrix  requires  only  a 
single  computer  word.  Besides  reducing  the  storage  require¬ 
ment,  packing  speeds  up  row  operations  since  logical  opera¬ 
tions  between  corresponding  bits  can  be  performed  simultaneously. 

Many  graphical  operations  can  be  reduced  to  a  suc¬ 
cession  of  row  operations.  For  example  a  fast  routine  for 
forming  the  transitive  closure  of  a  relation  represented  by 
a  Boolean  matrix  due  to  Warshall  [10]  is  the  following: 

Scan  the  Boolean  matrix  using  a  vertical  raster 
scan.  Whenever  a  1  is  encountered  in  row  i 
column  j,  replace  row  i  by  the  logical  inclu¬ 
sive  OR  function  of  row  i  and  row  j. 

This  routine  applied  to  the  matrix  of  Figure  19  results  in 

the  matrix  of  Figure  20.  The  matrix  SUC  represents  the 

transitive  closure  of  matrix  OUT. 

The  program  uses  the  matrix  SUC  to  calculate  the 
maximum  algorithm  subsets.  The  subroutine,  vrhich  forms 
maximum  algorithm  subsets,  first  calculates  the  condensation 
of  matrix  SUC  with  respect  to  matrix  ALG.  If  a  diagonal 
element  of  the  resulting  matrix  is  equal  to  1,  the  algorithm 
corresponding  to  the  diagonal  element  is  not  a  compatible 
set.  A  submatrix  of  SUC  containing  only  rows  and  columns 
corresponding  to  outputs  in  this  algorithm  is  formed.  The 
maximal  subsets  art  calculated  from  the  forward  and  converse 
level  partitions  of  the  subgraph  and  the  refined  maximal 


MATRIX  SOC 


1  2  3 

1 2345678901 2345678901 2345678901 

0000000000000000000000000000000 
0C001 01 01 0000000000001 1 00000000 

oooioi  oi  oi  ooi  i  on  oi  n  ooi  oi  on  oi 
0000000000000000000000000000000 
0000001010000000000001 1 00000000 
0001 000101 001 1 011 0111 001 01 01 1 01 
000000001 0000000000001 1 00000000 
oooi  oooooi  ooi  i  on  om  ooi  oi  oi  i  oi 
0000000000000000000001 1 00000000 
oooi  ooooooooi  i  on  on  i  ooi  oi  oi  i  oi 
1 000000000000000000000000000000 
0000000000000000000000000000000 
0000000000000000000000000000001 
0000000000001 000000000000000001 
oooi  oooi  oi  ooi  i  on  om  ooi  oi  oi  i  oi 
0001 000000001 1 00000000000001 001 
0001 000000001 1 01 000000000001 001 
0001 000001 001 10110111 001 01 01 1 01 
0001 000000001 1 01 1 0001 001 01 01 1 01 
0001 000000001 101101 01 001 01 01 1 01 
0001 000000001 1 01 1 000000001 01 1 01 
00000000000000000000001 00000000 
0000000000000000000000000000000 
0001 000000001 1 01 1 OOOI 00001 01 1 01 
000000000000000000000000001 0000 
0001 000000001 1 01 1 00000000001 1 01 
0000000000000000000000000000000 
0000000000001 1 00000000000000001 
0001 000000001 1 01 1 00000000001 001 
000000000001 0000000000000000000 
0000000000000000000000000000000 


FIGURE  20 

otal  Output  Dependency  Relation  Matrj 


subsets  are  substituted  for  the  original  algorithm  definition. 
The  resulting  maximum  subset  definitions  are  shown  in 
Figure  21.  These  splits  are  displayed  in  Figure  22.  The 
reader  will  note  that  the  subroutine  was  unable  to  refine  the 
split  of  algorithm  8.  Output  22  is  in  maximum  algorithm 
number  8  and  maximum  algorithm  number  18. 

The  condensation  of  the  network  relation  with 
respect  to  the  maximal  subsets  is  represented  by  the  Boolean 
matrix  of  Figure  23.  This  matrix  represents  the  maximal 
subset  output  algorithm  graph.  The  maximal  connected  sub¬ 
sets  of  this  graph  are  constructed  by  the  following  routine: 

1.  Form  the  transitive  closure  of  the  algorithm 
dependency  relation. 

2.  Obtain  the  symmetric  portion  of  this  matrix. 

3.  Set  diagonal  elements  equal  to  1. 

Each  distinct  row  of  the  resulting  matrix  is  a  maximal 
connected  subset.  The  maximal  connected  algorithm  subsets 
for  the  example  problem  are  shown  in  Figure  24  and  the 
dependency  relation  among  maximal  connected  algorithm  sub¬ 
sets  is  represented  by  the  Boolean  matrix  of  Figure  25. 

The  converse  level  partition  of  the  subset  dependency  re¬ 
lation  of  Figure  25  is  shown  in  Figure  26.  According  to 
this  partition  maximal  connected  algorithm  subset  11  occurs 
first  in  the  simulation  sequence.  This  subset  includes  only 
maximal  algorithm  13  which  is  a  subset  of  algorithm  1.  The 
first  algorithm  in  the  simulation  sequence  is  thus  algorithm  1 
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The  second  and  third  connected  algorithm  subsets  are  2  and 
13.  These  subsets  correspond  to  algorithms  2  and  6.  The 
subsets  3  and  15  correspond  to  algorithms  3  and  7*  Sub¬ 
sets  4,  16  and  14  correspond  to  algorithms  4,  8  and  7 
respectively.  Subset  8  contains  maximal  algorithms  8  and  12 
Maximal  algorithm  8  contains  outputs  21  and  22  while  maximal 
algorithm  12  contains  outputs  23,  24  and  25.  Output  22  has 
'’l"  simulation  state  since  it  was  updated  by  a  prior  applica 
tion  of  algorithm  8.  With  this  initial  condition  the  modi¬ 
fied  eligible  set  method  forms  the  simulation  sequence  12 
8  for  subset  8. 

Except  for  subset  5  the  remaining  subsets  corres¬ 
pond  to  single  algorithms.  The  modified  eligible  set 
method  forms  the  sequence  595  for  subset  5.  The  resulting 
simulation  sequence  is 

12637487  12  8  11  10  76  10  65951. 

The  only  searching  required  to  obtain  the  above 
solution  was  in  the  construction  of  a  simulation  sequence 
for  subset  5  and  only  two  sequences  were  examined  in  the 
construction  of  a  simulation  sequence  for  this  subset. 

The  modified  eligible  set  method  was  used  to 
calculate  a  simulation  sequence  for  this  network  without 
the  aid  of  algorithm  splitting.  This  method  used  a  maximum 
of  12  levels.  The  number  of  choices  at  each  level  varied 
from  2  to  6.  Approximately  500,000  potential  sequences  were 
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examined.  With  the  implementation  of  the  algorithm  splitting 
and  partitioning  techniques  the  problem  became  trivial. 
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VIII.  EXTENSIONS 

In  the  previous  chapters  a  sir.:;  le  characterization 
of  the  simulation  sequence  ordering  problem  was  employed. 

In  this  chapter  various  extensions  of  the  problem  are 
discussed . 

In  the  system  model,  it  was  assumed  that  each 
process  contained  a  single  algorithm  which  calculated  new 
values  for  the  set  of  state  variables.  This  assumption 
greatly  simplified  the  construction  of  simulation  sequences. 

In  Section  8.1  the  single  algorithm  per  set  of  state  vari¬ 
ables  assumption  is  replaced  by  a  single  algorithm  per 
process  assumption. 

In  the  previous  chapters  it  was  assumed  that  a 
minimum-length  simulation  sequence  was  desired.  Actually 
the  simulation  sequence  with  minimum  cost  is  usually  desired. 
In  Section  8.2  each  algorithm  is  assumed  to  have  a  positive 
weight  associated  with  it  and  the  formation  of  a  simulation 
sequence  with  minimum  linear  sum  is  discussed. 

In  Section  8.3,  the  modelling  of  random  variables 
is  discussed.  The  problem  of  forming  a  general  simulation 
sequence  for  a  network  with  several  configurations  is  dis¬ 
cussed  in  Section  8.4.  In  Section  8.5  the  elimination  of 
time  steps  is  discussed. 

8.1  Single  Algorithm  Characterization 

In  Chapter  II  in  order  to  guarantee  that  a  simula¬ 
tion  sequence  will  exist  for  any  process  it  was  assumed  that 
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a  single  algorithm  calculated  new  values  for  the  set  of 
state  variables.  In  this  section  it  will  be  assumed  that 
a  single  algorithm  exists  for  each  process.  This  algorithm 
will  calculate  new  values  for  all  process  output  and  state 
variables.  This  is  the  characterization  employed  by 
Parnas  [10].  Parnas  considers  two  types  of  application  of  a 
process  algorithm,  Hp.  In  the  first  type  of  application, 
denoted  by  Hp,  new  values  are  obtained  for  all  process  p 
output  and  state  variables.  In  the  second  type  of  applica¬ 
tion,  denoted  by  <vHp,  new  values  are  obtained  for  all  pro¬ 
cess  p  output  and  state  variables,  but  the  state  variable 
values  are  restored  to  their  initial  value  after  the  applica¬ 
tion  of  the  algorithm.  An  application  of  MIP  thus  effect¬ 
ively  calculates  new  values  only  for  output  valuables. 

Parnas  assumes  that  only  the  input-output  process 
dependency  relation  is  known.  Parnas  [10]  states:  "Correct 
values  for  the  state  variables  can  be  determined  if  and  only 
if  the  event  description  algorithm  is  executed  with  all  in¬ 
put  variables  known  to  be  correct."  He  therefore  assumes 
that  each  process  state  variable  is  dependent  upon  all  the 
inj:  t.  and  state  variables  of  this  process. 

Using  this  characterization  a  minimum-length 

simulation  sequence  can  be  obtained  as  follows: 

«  * 

1.  Add  to  each  process  an  output  y  .  The  value  of  y 
is  assumed  to  be  dependent  upon  all  process  inputs. 
Output  y*  is  not  related. 
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2.  Assume  tnat  for  each  process  there  exists  a  single 
algorithm  which  will  calculate  new  values  for  each 
process  output  (including  y*). 

3.  Obtain  a  minimum-length  output  simulation  sequence. 

4.  For  each  process  p,  scan  the  output  simulation 
sequence  for  occurrences  of  the  process  p  output 
algorithm.  Label  the  last  occurrence  HP  and  all 
preceding  occurrences  MlP. 

The  resulting  sequence  is  a  minimum-length  simulation  sequence 

for  the  single  algorithm  characterisation. 

Consider  the  example  problem  discussed  in  Section  2.3. 

The  network.  for  this  example  is  displayed  in  Figure  J|.  The 

network  output  graph  is  displayed  in  Figure  27.  A  minimum- 

length  output  simulation  sequence  is  F  FJF  F  F  .  A  minimum- 

length  simulation  sequence  for  the  single  algorithm  charac- 

2  3  12"? 

terization  is  thus  'V>H  H  H  H'  . 

8.2  Weighted  Algorithms 

In  the  previous  chapters  it  has  been  assumed  that 
a  minimum-length  simulation  sequence  is  desired.  If  the 
cost  of  each  algorithm  is  known,  a  simulation  sequence  with 
minimum  cost  is  usually  desired.  The  problem  of  constructing 
such  a  sequence  is  discussed  in  this  section.  Each  algorithm 
is  assuirec.  to  have  a  positive  weight  associated  witti  it  and 
the  simulation  sequence  with  minimum  linear  sum  is  desired. 

The  approach  employed  in  the  formation  of  a 
minimum-length  simulation  sequence  is  to  include  each  algo¬ 
rithm  a  minimum  number  of  times  in  the  simulation  sequence. 
Since  each  algorithm  has  positive  weight,  such  a  sequence 
is  a  minimum-weight  simulation  sequence  for  any  positive 


FIGURE  27 

MODIFIED  NETWORK  OUTPUT  GRAPH 
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weight  assignment.  The  algorithm  splitting  technique  dis¬ 
cussed  in  Chapter  VI  is  thus  applicable  to  this  problem. 

In  Chapter  V  a  forward  level  partition  of  the 
algorithm  graph  was  shown  to  form  a  minimum-length  simulation 
sequence.  In  this  sequence  each  algorithm  is  included  exactly 
once.  This  sequence,  is  also  a  minimum-weight  sequence  for 
any  positive  weight  assignment.  Furthermore,  in  the  case  of 
maximal-connected  subsets,  an  ordering  of  the  minimum-weignt 
simulation  sequences  for  the  maximal- connected  subsets  will 
form  a  minimum-weight  simulation  sequence. 

The  only  change  required  is  In  the  algorithm  of 
Section  7.2.  This  algorithm  uses  a  selective  search  techni¬ 
que.  The  reader  will  note  that  the  optimal  choice  criteria 
is  still  applicable.  The  potential  length  test  must,  how¬ 
ever,  be  replaced  by  a  potential  weighted  sum  test.  The 
weighted  algorithm  characterization  thus  requires  only  a 
trivial  modification. 

8.3  Random  Variables 

Random  variables  are  frequently  employed  in  simu¬ 
lation  problems.  In  the  system  model  it  was  assumed  that 
each  function  is  deterministic.  In  this  section  the  modelling 
of  random  variables  is  discussed. 

Associated  with  each  random  variable  is  a  function 
h  such  that  successive  applications  of  the  function  h  will 
yield  a  set  of  random  values  In  the  desired  distribution. 

Each  application  of  the  function  will  result  in  a  different 


value. 
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The  random  variable  is  assumed  to  be  distributed 
over  simulated  time.  It,  therefore,  has  only  a  single  value 
associated  with  each  time  step,  and  all  variables  depend¬ 
ing  upon  a  random  variable  should  use  the  same  value  for 
the  random  variable.  This  can  be  achieved  by  considering 
each  random  variable  to  be  a  state  variable  that  is  R-re- 
lated  to  itself.  If  the  random  variable  supplies  values  to 
variables  in  different  processes,  an  output  variable  whose 
value  is  equal  to  the  random  variable  can  be  added  to  the 

process.  By  the  definition  of  a  simulation  sequence  a 
function  associated  with  a  self-related  state  variable 
must  be  applied  exactly  once. 

8.4  Time-Varying  Network  Structure 

Suppose  that  the  network  structure  of  a  simulation 
problem  assumes  two  or  more  configurations  for  various  time 
steps  of  the  simulation  and  that  a  single  simulation  sequence 
is  desired.  This  sequence  must  be  a  simulation  sequence  for 
each  configuration  of  the  network .  In  this  section  a  techni¬ 
que  for  constructing  a  general  simulation  sequence  for  a  multi 
configuration  network  is  presented. 

Each  configuration  of  the  network  is  assumed  to 
be  represented  by  an  acyclic  network  graph.  These  graphs 
are  assumed  to  be  disjoint  and  subgraphs  of  a  general  network 
output  graph.  The  general  network  output  graph  contains  a 
set  of  nodes  for  each  output.  The  application  of  an 
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algorithm  containing  an  output  will  calculate  new  values  for 
each  output  in  the  set  of  nodes  corresponding  to  this  output. 

A  simulation  sequence  for  this  general  network  output  graph 
is  a  simulation  sequence  for  each  configuration  of  the 
network , 

The  network  output  grains  for  three  configurations 
of  a  network  are  shown  in  Figure  28.  The  general  network 
output  graph  for  this  example  is  shown  in  Figure  29.  The 

i  p  -3  I 

sequence  F  F  FJF  is  a  minimum-length  simulation  sequence 

for  the  general  network  output  graph.  This  sequence  is  a 

simulation  sequence  for  each  configuration  of  the  network. 

A  minimum-length  simulation  sequence  for  a  general 

network  output  graph  is  the  shortest  seque  ls  which  is  a 

simulation  sequence  for  every  configuration  of  the  network. 

Some  configurations  may  have  a  shorter  simulation  sequence. 

For  the  configuration  containing  outputs  21  through  26, the 
l  p  3 

sequence  F  F  FJ  is  a  simulation  sequence. 

The  advantage  of  a  general  simulation  sequence 
for  a  set  of  network  configurations  is  that  separate  simula¬ 
tion  sequences  need  not  be  stored  for  each  configuration  and  that 
a  selection  of  sequence  need  not  be  performed.  The  disad¬ 
vantage  of  a  general  simulation  sequence  is  that  this 
sequence  is  in  general  not  a  minimum-length  sequence  for 
every  configuration. 


FIGURE  28 


NETWORK  OUTPUT  GRAPHS  FOR 
THREE  NETWORK  CONFIGURATIONS 
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FIGURE  29 


GENERAL  NETWORK  OUTPUT  GRAPH 


8.5  Time-Step  Minimization 

The  problem  of  constructing  a  simulation  sequence 
for  a  time-step  has  been  discussed.  A  simulation  problem 
generally  has  a  large  number  of  time-steps.  In  this  section 
the  reduction  of  the  number  of  time-steps  is  discussed. 

Suppose  that  each  output  has  a  small  amount  of 
delay  associated  with  it  and  that  the  network  is  quiescent . 
Due  to  the  small  amount  of  delay  associated  with  each  output, 
any  change  of  value  in  an  input  or  state  variable,  called 
an  event ,  cannot  propagate  instantaneously  through  the  net¬ 
work  but  will  ripple  through  the  network.  If  the  network 
is  acyclic  and  no  new  event  occurs,  the  outputs  will  reach 
a  quiescent  state. 

In  general  many  time  steps  are  required  before 
the  quiescent  state  is  reached.  If  the  transient  behavior 
of  the  network  is  not  required,  these  time  steps  can  be 
eliminated  by  calculating  the  quiescent  state  output  values 
directly. 

A  quiescent  value  can  be  calculated  for  an  output 
only  if  the  outputs  on  which  this  output  depends  have  quie¬ 
scent  value.  This  dependency  relation  is  the  same  relation 
as  the  output  dependency  relation  used  In  the  formation  of 
a  simulation  sequence.  The  output  simulation  sequence 
will  thus  calculate  the  quiescent  state  values. 
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The  degree  to  which  the  reduction  of  time  steps 
can  be  employed  is  dependent  upon  the  output  desired  from 
the  simulation  and  the  occurrence  of  events  within  the 
simulation. 
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IX.  SUMMARY  AND  CONCLUSIONS 

9.1  General 

The  problem  of  forming  a  simulation  sequence  for 
networks  with  interactive  processes  has  been  discussed.  The 
process  interaction  is  specified  in  the  system  model  by  the 
process  and  network  relations.  The  problem  of  forming  a  simu¬ 
lation  sequence  for  the  system  model  was  reduced  to  the  problem 
of  forming  an  output  simulation  sequence. 

An  output  dependency  relation  was  defined  and  repre¬ 
sented  by  the  network  output  graph.  A  simulation  sequence 
exists  only  for  systems  represented  by  an  acyclic  network 
output  graph.  Paths  in  the  network  output  graph  correspond  to 
subsequences  of  simulation  sequences.  Partial  orderings  of  the 
set  of  outputs  were  represented  by  the  forward  and  backward 
level  sets. 

An  algorithm-dependency  relation  was  obtained  from 
the  output  dependency  relation  via  a  condensation  with  respect 
to  the  algorithm  definition.  This  relation  is  represented  by 
the  algorithm  graph.  If  the  relation  is  acyclic,  it  specifies 
a  minimum-length  simulation  sequence.  If  the  relation  is 
cyclic,  an  ordering  of  the  set  of  algorithms  is  not  possible. 

An  ordering  of  the  maximal  connected  subsets  of  the  algorithm 
graph  is,  however,  always  possible.  Furthermore,  a  concatena- 
-  tion  of  minimum-length  simulation  sequences  for  maximal  con¬ 
nected  subsets  was  shown  to  form  a  minimum-length  simulation 
sequence. 
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A  criteria  for  splitting  algorithms  has  been  presented. 
Algorithm  which  are  not  compatible  sets  must  be  multiple  in¬ 
cluded  in  simulation  sequences  and  hence  can  be  partitioned. 

A  general  level  partition  of  an  algorithm  into  maximal  subsets 
and  refinements  of  these  subsets  has  been  described. 

A  selective  search  technique  for  the  construction  of 
simulation  sequences  for  problems  represented  by  a  connected 
algorithm  graph  of  maximal  algorithm  subsets  has  been  presented. 
This  algorithm  employs  a  one-step  look  ahead,  recognizes  optimal 
choices  and  includes  a  potential  length  test. 

A  general  method  for  the  construction  of  simulation 
sequences  was  described.  This  method  uses  all  the  techniques 
developed. 

By  means  of  simple  modifications  this  method  can  be 
used  to  form  simulation  sequences  for  the  single  algorithm  per 
process  characterization  of  the  problem.  The  construction  of 
minimum-cost  simulation  sequences  for  a  network  with  positive 
algorithm  costs  required  only  a  change  in  the  potential  length 
test  employed  by  the  selective  search  routine. 

9.2  Contributions  of  the  Research 

The  contributions  of  this  research  are  mainly  in  the 
area  of  digital  computer  simulations  and  especially  relating  to 
the  simulation  of  systems  composed  of  multiple-output  interactive 
subsystems.  The  following  list  is  a  summary  of  the  significant 


contributions : 


1.  The  formation  of  simulation  sequences  has  been  shown 
to  be  independent  of  process  state  variables.  Only 
the  output  dependency  relation  need  be  considered. 

2.  The  algorithm  dependency  relation  has  been  shown  to 
partition  the  simulation  sequence  formation  problem 
Into  a  collection  of  smaller  problems. 

3.  An  algorithm  splitting  technique  has  been  described. 
This  technique  greatly  reduces  the  amount  of  compu¬ 
tation  in  the  formation  of  simulation  sequences  for 
many  problems. 

4.  A  general  method  for  the  formation  of  minimum-length 
simulation  sequences  employing  (1),  (2),  (3)  and  a 
selective  search  routine  has  been  presented. 

9 . 3  Future  Extensions 

The  work  described  in  this  thesis  has  suggested  the 
following  areas  for  further  research. 

1.  Development  of  an  analytic  routine  for  forming  simula¬ 
tion  sequences  for  a  connected  algorithm  graph  composed 
of  maximal  algorithm  subsets. 

2.  Investigation  into  the  interaction  of  processes  over 
time  intervals  greater  than  a  single  time-step. 

3.  Examination  of  algorithm  equivalences  and  the  examina¬ 
tion  of  simulation  sequence  equivalences. 
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Prepared  by:  Dieter  Knollman 
Applicable  to:  Univac  1108 
Type :  Software 


New  York  University 
Department  of  Electrical  Engineering 
Computer  Bulletin 

SIMULATION  SEQUENCE  FORMATION  PROGRAM 

I.  ABSTRACT 

This  program  will  form  a  minimum-length  simulation 
sequence  for  an  interactive  discrete-time  system.  The  system 
is  represented  by  a  block  diagram.  The  input  to  the  program 
consists  of  a  specification  of  the  output  dependency  relation 
and  a  listing  of  outputs  included  in  each  algorithm.  The 
output  of  the  program  is  a  minimum-length  string  of  algorithms 
that  will  accurately  simulate  the  behavior  of  the  system. 

II.  PERFORMANCE  SPECIFICATION 

The  program  can  not  form  simulation  sequences  for 
problems  with  cyclic  output  dependencies.  For  these  problems 
a  listing  of  all  cyclic  output  dependencies  is  produced. 

The  program  requires  that  all  outputs  and  algorithms 
are  numbered  in  consecutive  order  using  one-indexing.  Empty 
algorithms  are  allowed.  Each  output  must  be  included  in  a 
single  algorithm.  If  this  requirement  is  violated,  the 
program  will  produce  a  list  of  outputs  that  are  included 
in  no  algorithm  and  a  list  of  outputs  that  are  included  in  two 
or  more  algorithms. 
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The  speed  of  the  program  is  problem  dependent. 

The  core  requirements  are  determined  by  the  number  of  out¬ 
puts.  The  parameter  IPA  determines  the  size  of  data  arrays. 
This  parameter  must  be  greater  than  the  number  of  outputs. 

In  the  program  IPA  is  equal  to  5Q.  IPA  can  be  altered  by 
changing  a  single  parameter  definition  statement. 

III.  DESCRIPTION 

The  program  consists  of  a  collection  of  external 
subroutines  and  a  main  program.  The  main  program  determines 
the  storage  allocation  and  schedules  the  execution  of  sub¬ 
routines.  A  description  of  the  main  program  and  the  sub¬ 
routines  is  presented  in  the  listing. 

IV.  OPERATING  INSTRUCTIONS 

The  necessary  control  cards  to  compile  and  execute 
the  program  are  supplied  with  the  program.  The  user  must 
supply  a  RUN  card  and  data  cards.  The  algorithm  definition 
cards  must  precede  the  output  dependency  definition  cards. 
The  format  of  the  data  cards  is  presented  in  the  listing  of 
the  READA  and  the  READK  subroutines . 

V.  EXAMPLE 

A  typical  run  will  consist  of  the  following: 

RUN  Card 

FORTRAN  V  deck  with  control  cards 
Algorithm  definition  data  cards 
Output  dependency  data  cards 
EOF  card 
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